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Preface 

AFCRL  support  for  the  SAMS/ABRES  Rain-Erosion  Project  at  WallopB 
Island,  Virginia,  during  the  1970-73  period  covered  in  this  report,  was  furnished 
primarily  by  Air  Force  and  contractor  personnel  of  the  Meteorology  Laboratory 
(LY),  Directed  by  Dr.  Morton  L.  Barad.  The  Laboratory  Manager  for  Weather 
Erosion  Programs  was  Chankey  N.  Touart  (from  June  1973  to  present).  The 
Wallops  Island  Measurement  Program  was  directed  by  Dr.  Robert  M.  Cunningham. 
Tne  radar  measurement  aspects  of  the  program  were  supervised  by  Dr.  Kenneth 
R.  Hardy. 

The  AFCRL  and  contract  contributors  to  the  radar  measurement  program 
were  identified  in  AFCRL/SAMS  Report  No.  1.  These  persons  also  contributed, 
directly  or  indirectly,  to  the  work  reported  herein. 

Personnel  associated  with  the  aircraft  measurement  program  supporting 
SAMS/ABRES  at  Wallops  Island  have  not  been  identified  previously.  The  follow¬ 
ing  list  of  personnel  and  their  work  efforts  toward  the  accomplishments  of  the 
objectives  of  the  rain-erosion  program  are  acknowledged.  Particular  thanks  are 
also  expressed  to  Hugh  Church  of  the  Sandia  Corporation,  Albuquerque,  N.  M. , 
and  to  Alfred  A.  Spatola,  I.t  Col  James  F.  Church,  and  Elizabeth  L.  Kintigh  of 
AFCRL  for  their  reviews  and  criticisms  of  the  manuscript  of  this  report.  Their 
suggestions  were  most  helpful  and  appreciated. 
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*I.YC  Convective  Cloud  Physics  Branch,  Meteorology  Laboratory,  AFCRL. 
**MRI  Meteorology  Research  Inc. ,  Altadena,  California. 

The  SAMS  (Sandia  Air  Force  Materials  Study )  Program  began  in  1069  as  a 
jointly  funded  effort  between  SAMSO  'Air  Fnm*  and  Sandia  I  .aboratories.  By  1973 
all  funding  was  being  done  bv  SAMSO.  The  objective  of  the  SAMS  Program  is  to 
experimentally  test  the  effects  df  precipitation  and  cloud  particle  (hydrometeor) 
impacts  upon  various  full  scale  missile  materials  by  Hying  high  speed  vehicles 
through  actual  storm  and  cloud  environments.  The  test  vehicle  is  launched  at  a 
relatively  low  elevation  angle  (typically  30°)  and  performs  the  impact  erosion 
experiment  on  the  ascending  portion  of  its  trajectory,  The  instrumented  payload 
with  its  test  nnsetip  and  heatshield  is  subsequently  recovered  from  the  ocean  by 
means  of  a  parachute  and  notation  system. 

The  NASA  Wallops  Flight  Center  on  the  eastern  shore  of  Virginia  was  selected 
for  the  site  of  these  tests  because  it  exhibits  a  relathelv  high  frequency  of 
occurrence  of  widespread  stratiform  storms,  has  the  necessary  support  facilities, 
and  is  readily  accessible.  Storm  Invironinents  are  measured  1>\  instrumented 
aircraft  and  ground  inst  i  i merit s  ineluding  special  weather  radars.  Further  details 
of  the  test  set  up  a  in  .  outlined  in  the  SA  \IS  Program  Test  Plan.  I>\  J.K.  Cole, 
S('-l)R-70-Hr>0,  Sand  in  1  ,n  born  to  ties .  Mbuqiierque.  New  Mexico  December  1970, 
while  test  results  are  contained  in  \nrious  i  iissifieil  reports  from  Sandia 
I  .allocator  ies. 
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Hydrometeor  Parameters  Determined  from  the 
Radar  Data  of  the  Sams  Rain  Erosion  Program - 

AFCRL/Sams  Report  No.  2 


1.  INTRODUCTION 

The  radar  equations  and  measurement  techniques  of  the  SAMS  rain-erosion 

program  at  Wallops  Flight  Center,  Virginia  have  been  described  in  a  previous  re- 
1  * 

port  by  Plank.  ’  It  was  noted  therein  that  the  end-product  radar  data  for  the 
SAMS  missile  trajectories  consist  of  measurements  of  the  radar  reflectivity  fac¬ 
tor  for  water  hydrometeors,  Zw ,  in  the  lowest  storm  portion,  containing  rain, 
and  of  the  radar  reflectivity  factor  for  ice  hydrometeors,  Zj,  in  the  upper  storm 
regions,  containing  snow  and  ice  crystals.  It  was  also  mentioned  that  the  values 
of  the  reflectivity  factor  are  indeterminant  within  the  melting  zone  in  the  storm, 
since  the  radar  backscattering  properties  of  an  ice-water  mixture  of  falling, 
melting  snowflakes  are  not  describable  in  general  quantitative  terms,  at  least,  not 
at  present. 

The  purpose  of  this  report  is  to  explain  how  the  radar  values  of  Z^y  and  Zj 
for  the  SAMS  missile  trajectories  were  used  to  ascertain  the  corresponding 
values  of  the  precipitation  rate,  P,  the  precipitation  liquid- water- content,  M, 
(Received  for  publication  4  June  1974) 

1.  Plank,  V.G.  (1974)  A  summary  of  the  radar  equations  and  measurement 
techniques  used  in  the  SAMS  rain  erosion  program  at  Wallops  Island, 
Virginia,  AFCRL/SAMS  Report  No,  1,  108  pps,  AFCRL-TR-74-0053, 

Special  Reports  No.  172,  Air  F*orce  Cambridge  Research  Laboratories, 
Bedford,  Mass. 

-The  name  "Wallops  Flight  Center"  was  officially  designated  in  May  1974.  The 
facility,  prior  to  this,  was  called  "Wallops  Station  . 
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and  other  related  parameters  of  SAMS  interest,  such  as  the  size-distribution 
properties  of  the  hydrometeors. 


2.  BACKGROUND  OF  REFLECTIVITY- FACTOR 
VS  PRECIPITATION-RATE  RELATIONSHIPS 

2.1  Relation*  for  Rain 

2  3  4 

W'exler,  Marshall  et  al,'  and  Marshall  and  Palmer  were  the  first  to  inves¬ 
tigate  the  theoretical-experimental  relationships  between  the  back- scattered  sig¬ 
nals  detected  by  radar  and  the  hydrometeor  parameters  cited  above.  Their  spe¬ 
cific  interest  involved  the  possibility  of  using  radar  to  determine  the  intensity  and 

distribution  of  rainfall  for  hydrological  purposes.  Wexler,  from  the  data  of  Laws 
5 

and  Parsons,  and  Marshall  et  al,  from  their  own  acquired  data,  ascertained  that 
the  radar  reflectivity  factor  for  water,  Z^. ,  was  related  to  the  precipitation  rate, 
P ,  approximately  as  described  by  the  equation 

Zw  =  kpe.  (1) 

”1  6  —  3 

where  P  is  specified  (usually!  in  mm  hr  and  Z w  is  given  in  mm  m 

Wexler  found  that  K  had  the  value  320,  with  E  =  1.44.  Marshall  et  al^  deter- 

4 

mined  that  K  =  190  and  E  =  1.72.  Marshall  and  Palmer'  presented  revised  val¬ 
ues  for  the  Marshall  et  al  data  of  K  =  220  and  E  =  1.60.  They  also  reported  val¬ 
ues,  based  on  other  data,  of  K  =  296  and  E  =  1.47  . 

Many  subsequent  investigators  have  reported  particular  K  and  E  values  for 

various  rainfall  situations  of  experimental  observation.  These  have  been  reviev'ed 

6  7  8  9 

and  listed  by  Sivaramakrishnan,  Atlas,  Stout  and  Mueller,  Battan,'  and  others. 


2.  Wexler,  R.  (1947)  Radar  detection  of  a  frontal  storm,  18  .June  1946, 

Meteorol.  4:38-44. 

■  -  -  ■  «*v 

3.  Marshall,  J.  S.  ,  I.angille,  R.C.,  and  Palmer,  W\  McK.  (1947)  Measurement 

of  rainfall  by  radar,  J.  Meteorol.  4(6):  186-192, 

i  ■  —  —  —  vw 

4.  Marshall,  J.S. ,  and  Palmer,  W.  McK.  (1948)  The  distribution  of  raindrops 

with  size,  .1,  Meteorol.  5:165-166. 

—  —  ■  —  VMr 

5.  Laws,  .1.0.,  and  Parsons,  D.A.  (1943)  The  relation  of  raindrop  size  to  in¬ 

tensity,  Trans.  Amer.  Geophys,  Union  24(Part  II):452-460. 

6.  Sivaramakrishnan,  M.V.  (1961)  Studies  of  raindrop  size  characteristics  in 

different  types  of  tropical  rain  using  a  simple  raindrop  recorder,  Indian 
J.  Met,  Geophys.  12G89  (473-9,  602,  609,  613). 

7.  Atlas,  D.  (1964)  Advances  in  radar  meteorology,  Adv. Geophys.  10:317  (399, 

434,441,449,478,481). 

8.  Stout,  G.  E. ,  and  Mueller,  E.  A.  (1968)  Survey  of  relationships  between  rain¬ 

fall  rate  and  radar  reflectivity  in  the  measurement  of  precipitation,  J. 
Appl.  Meteorol.  L465  (478-9).' 

9.  Battan,  L..I.  (1973)  Radar  Observation  of  the  Atmosphere,  the  University  of 

Chicago  Press,  Chicago,  Ill. 
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Battan,  for  example,  lists  69  sets  of  values  for  different  investigators,  geo¬ 
graphical  locations,  seasons,  and  rainfall  type  (widespread,  orographic,  showery). 
The  K  values  range  from  16.6  to  730;  the  E  values  from  1.18  to  2.87. 

2.2  Relation!  for  Snow 

Marshall  and  Gunn,1®  fro-"  the  data  of  Langille  and  Thain,11  found  that  a 
relationship  of  the  Eq.(l)  type  could  also  be  applied  to  the  case  of  snow  crystals 
falling  in  the  atmosphere.  Measurements  of  Zj  were  obtained  by  radar  and  were 
correlated  with  values  of  P ,  computed  from  snowflake  samples  collected  on  filter 
paper  during  the  same  time  period  as  the  radar  measurements.  They  found  that 

4 

the  Marshall  and  Palmer  equation  for  rain,  with  K  =  220  and  E  =  1.6,  provided 

a  good  description  of  the  data.  The  data  were  "best  fitted"  by  an  equation  in  which 

1 9 

K  =  200  and  E  =  1.6  .  Imai  et  al  likewise  found  that  a  relationship  of  the  Eq.  (1) 
type  was  descriptive  of  small  snow  crystals;  their  best  fit  values  were  K  =  500 
and  E  =  1.6.  Gunn  and  Marshall,1^' 1(1  reported  that  K  =  2000,  with  E  =  2.0  for 
aggregated  snow.  They  additionally  mentioned  that  the  snow  equation  presented 
earlier  by  Marshall  and  Gunn, 10  was  probably  erroneous,  because  of  problems 

15 

with  the  sampling  methods  used  at  that  time.  [However,  Carlson  and  Marshall, 

with  reference  to  private  communication  by  K.  L.  S.  Gunn,  have  indicated  that  the 

1 6 

Marshall  and  Gunn  equation  is  descriptive  of  single  crystal  snow.]  Austin 

cited  values  of  K  =  1000  and  E  =  1.6,  for  heavy  snow  in  the  New  England  area. 

Her  values  were  intermediate  between  those  of  Marshall  and  Gunn10  and  those  of 
13  14 

Gunn  and  Marshall. 


10.  Marshall,  J.  S. ,  and  Gunn,  K.  L.  3.  (1952)  Measurement  of  snow  parameters 

by  radar,  J.  Meteorol.  9:322. 

11.  Langille,  R.  C. ,  and  Thain,  R.  S.  (1951)  Some  quantitative  measurements  of 

three  centimeter  radar  echoes  from  falling  snow,  Canadian  J.  Phys.  29:482. 

12.  Imai,  I.,  Fujiwara,  M. ,  Ichimura,  I.,  and  Toyama,  Y.  (1955)  Radar  re¬ 

flectivity  of  falling  snow.  Pap,  in  Meteorol.  and  Geophys.  (.Tapan)J>:  130- 139. 

13.  Gunn,  K.  L.  S. ,  and  Marshall,  . I.  S.  (1956)  Size  distributions  of  aggregate 

snowflakes,  Scientific  Report  NW-20/B,  USAF  Contract  AF19(122)-217, 
McGill  Univ. ,  Montreal,  Canada. 

14.  Gunn,  K.  L.  S. ,  and  Marshall,  .  1.  S.  (1958)  The  distribution  with  size  of  ag¬ 

gregate  snowflakes,  J.  Meteorol.  1^5:452  (479). 

15.  Carlson,  P.  E. ,  and  Marshall,  J.  S.  (1972)  Measurement  of  snowfall  by  radar, 

J.  Appl,  Meteorol.  U:494-500. 

16.  Austin,  P.  M.  (1953)  Radar  measurements  of  the  distribution  of  precipitation 

in  New  England  storms,  Proc.  10th  Weather  Radar  Conf. ,  Washington, 
p.  247  (479.  491-482). 
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Ohtake  and  Henmi  determined  the  Zj  vs  P  relationships  for  various  types 
of  aggregated  snowflakes,  nine  in  all,  using  existing  .Japanese  data  concerning 
size  distribution,  particle  density,  and  particle  fall  speeds.  Their  results,  sum¬ 
marized  in  Table  l,  give  K  values  ranging  from  330  to  3300  and  E  values  rang¬ 
ing  from  1.5  to  2.3 . 


Table  1.  Examples  of  Zj  Vs  P  Relationships  for 
Different  Types  of  Aggregated  Snowflakes,  as  Re¬ 
ported  by  Ohtake  and  Henmi, Zj  is  in  mm**  m"9 
and  P  is  in  mm  hr"  * 


Type  of 

Aggregated 

Snowflakes 

Z  Vs  P 
Relation 

1.  Spatial  dendrites 

Zj  =  3300  P1,7 

2.  Plane  dendrites 

Zj  =  2900  P1-6 

3.  Stellars 

Zj  it  1800  P1-5 

4.  Thick  Plates 

Zj  =  1300  P2-3 

5.  Needles 

Zj  =  930  P1-9 

6,  Grauple 

Zj  =  900  P1-6 

7.  Bullets 

Zj  =  430  P1-5 

8.  Plates  and  Columns 

Zj  =  400  P1,6 

9.  Hail 

Z  =  330  P1,6 

i 

Additional  information  about  snow  of  different  types  has  been  presented  by 

18  19  20 

authors  such  as  I.angleben,  Kodairn  and  Inaba,  '  Litvinov,  Bashkirova  and 

21  22  23  15 

Pershina,  Kikuchi,  Sekhon  and  Srivastave,  and  Carlson  and  Marshall.1, 


17.  Ohtake,  T.  ,  and  Henmi,  T.  (1970'  Radar  Reflectivity  of  Aggregated  Snow¬ 

flakes,  preprints  of  papers  presented  at  the  1 4th  ftadar  Meteorology  Conf., 
Tucson,  Arizona,  17-20  November  1970,  pp.  209-211. 

18.  I.angleben,  M.  P.  (1954)  Theterminalvelocityofsnowaggrcgat.es,  Quart.  .1. 

Meteorol.  Soc.  80:174-181. 

,  - - - -  WV 

19.  Kodaira,  N.  ,  and  Inaba,  M.  (1955)  Measurement  of  snowfall  intensity  by 

radar.  Papers  .Meteorol, Geophys.  (.lapan)  J5:126-129. 

20.  Litvinov,  I.  V.  (1956)  Determination  of  the  terminal  velocity  of  snowflakes 

(in  Russian),  IZV.  Akad.  Nauk  SSSR  Ser,  Geofiz.  No.  7,  p.  853  (242). 

21.  Bashkirova,  G.  ,  and  Pershina,  T.  (1964)  On  the  mass  of  snow  crystals  and 

their  fall  velocities,  Tr,  GL  Geofiz,  Observ.  No.  165,  33-100. 

22.  Kikuchi,  K.  (1968)  On  snow  crystals  of  bullet  type,  .1.  Meteorol.  Hoc.  Japan  46: 

128-132.  '  ~~ 

23.  Sekhon,  R.S.  ,  and  Srivastava,  R.C.  (1970)  Snow  size  spectra  and  radar  re¬ 

flectivity,  .1.  Almos.  Sci.  27:209-307. 
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2.3  Relation*  lor  Ice  Cry  Main 

Information  about  the  Zj  vs  P  relationships  for  single- crystal  (non-aggre- 
gated)  ice  crystals  has  been  virtually  non-existent  in  the  literature  until  recently. 
This  is  primarily  (1)  because  ice  crystals,  such  as  those  of  SAMS  interest  that 
occur  in  the  uppermost  portions  of  storms,  seldom  exist  in  the  same  form  aloft 
as  at  the  ground  level,  where  most  previous  measurements  were  made,  and 
(2)  because  aircraft  instruments  suitable  for  measuring  the  ice  crystal  types  and 
properties  in  the  storms  themselves  have  not  been  developed  until  recently. 

Some  aircraft  measurement  data  have  been  acquired  since  1970  by  the  Uni¬ 
versity  of  Chicago  Cloud  Physics  Laboratory.  Optical  array  spectrometers  and 
formvar  replicator  instruments  were  used  to  determine  the  size,  type,  and  num¬ 
ber  concentration  of  ice  crystals  in  several  cirrus-cloud  and  contrail  situations, 

24  2*  26  27 

as  described  by  Kncllenberg  ’  ‘  *  and  Heymsfield  and  Knollenberg. 

Cunningham,  in  an  AFCRL  report  of  the  present  series  that  will  be  submitted 
later,  has  analyzed  two  sets  of  ice  crystal  data  supplied  to  AFCRL  by  Heymsfield. 
The  data  pertain  to  distribution  samples  of  ice  crystals  containing  columns,  bullets, 
and  rosettes.  Cunningham's  best  fit  curves  (of  Eq.  (1)  type)  for  the  computed  Zj 
and  P  values  yielded  K  =  16.1  and  E  =  1.69  for  the  first  data  set.  and  K  =  24.5, 
with  E  =  1.39  for  the  second. 

2.1  General  Application*  lo  SAMS 

Thus  far  in  this  background  review  the  Z ^  vs  P  relations  for  rain  and  the 
Zj  vs  P  relations  for  snow  end  ice  crystals  have  been  discussed.  Typical  values 
and  ranges  of  the  constants  and  exponents  of  Eq.  (1)  have  been  indicated. 

Empirical  equations  of  the  above  form  were  used  in  the  SAMS  rain-erosion 
program  to  determine  the  precipitation  rates  along  the  missile  trajectories  which 
corresponded  to  the  radar-measured  values  of  Z^  and  Zj .  Particular  hydro¬ 
meteor  regions  were  defined  for  the  Wallops  storm  experiments  (rain,  large-snow, 
small-snow,  and  ice-crystal  regions)  and  particular  values  of  the  coefficient,  K, 
and  exponent,  E,  of  Eq.(l)  were  used,  as  will  be  described.  Other  equations  of 


24.  Knollenberg,  R.  G.  (1970)  The  optical  array:  An  alternative  to  scattering  or 

extinction  for  airborne  particle  size  determination,  J.  Appl.  Meteor.  9:86- 
103. 

25.  Knollenberg,  R.G.  (1972)  Measurements  of  the  growth  of  the  ice  budget  in  a 

persisting  contrail,  .1,  Atmos.  Sci.  29:1367-1374. 

26.  Knollenberg,  R.G.  (1973)  Cirrus-contrail  cloud  spectra  studies  with  the 

Sabreliner,  Atmos.  Tech. ,  National  Center  for  Atmospheric  Research, 

No.  1,  March  1973,  pps.  52-55. 

27.  Heymsfield,  A. . I. ,  and  Knollenberg,  R.G.  (1972)  Properties  of  cirrus¬ 

generating  cells,  .T.  Atmos.  Sci.  29(7),  1358-1366. 
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power  function  form  were  also  used  to  relate  the  liquid- water-content  values,  M  , 
to  the  precipitation  rate,  and  to  interrelate  the  reflectivity  factors,  Zw  and  Z( , 
with  the  liquid-water  content.  These  latter  equations  will  be  presented  in  a  follow¬ 
ing  suction. 


3.  ADDITIONAL  IIACKCROIJND 

There  are  throe  topics  that  require  perspective  discussion  at  this  point  be¬ 
fore  the  specifics  of  the  AFCRL  analysis  program  for  the  SAMS  missile  shots  are 
described.  These  topics  concern: 

(H  the  fact  that  there  are  three  methods  of  evaluating  the  radar  reflectivity 
factor  that  require  identification,  one  defined  in  terms  of  radar  back-scattering 
theory,  one  defined  in  terms  of  the  actual  size-distribution  properties  of  the  hydro¬ 
meteors,  and  one  based  on  theoretical  distribution  functions  that  are  presumed  to 
be  descriptive  of  the  actual  distributions; 

(2)  the  definitions  of  precipitation  rate  and  liquid  water  content,  together  with 
discussion  of  how  these  parameters  are  measured  and  interrelated,  both  at  the 
surface  level  and  at  storm  levels  aloft, 

(3)  the  nature  .  >f  previous  correlation  and  r<  gression  analyses  and  comments 
about  sampling  volumes. 

These  topics  will  be  discussed  in  turn  in  the  following  sections. 


3.1  Different  Methods  of  Evaluating  the  Radar  Reflectivity  Factors 

As  noted  in  the  preceding  report  (Plank*),  a  weather  radar  basically  measures 
the  summed,  backscatter  cross-sections  of  the  hydrometeors  that  are  present 
within  the  radar  pulse  volume  (or  integration  volume,  if  video  integration  is  used). 
Specifically,  the  radar  measures  the  volume  reflectivity,  T? ,  which  is  the  summed 
backscatter  return  per  unit  volume.  The  radar  reflectivity  factors,  Z^, ,  for 
water  hydrometeors,  and  Zj ,  for  ice  hydrometeors,  are  then  determined  from 
the  equations  (Kqs.  (70)  and  (71)  of  the  cited  report). 


_  1.08  A  r?  v  ,  n12 


and 


z,  =  4-78rA  n  x  i°12  • 


(2) 


(3) 
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where  A  is  the  radar  wavelength,  in  cm ,  tj  is  specified  in  cm  * ,  and  and  Zj 
are  in  units  of  mm®  m*®  . 

It  is  important  to  note  that  no  assumptions  about  the  size  distribution  proper¬ 
ties  of  the  hydrometeors  are  made  when  Z^  and  Zj  are  computed  from  Eqs.(2) 
and  (3),  Such  assumptions  are  unnecessary,  since  the  summed  backscatter  re¬ 
turn  from  the  hydrometeors  is  being  measured  directly.  The  only  hydrometeor 
requirement  in  the  case  of  these  equations  is  that  information  be  available  to  dif¬ 
ferentiate  the  atmospheric  regions  containing  water  hydrometeors  from  those 
containing  ice  hydrometeors;  hence  to  establish  which  equation,  (2)  or  (3),  applies. 

For  future  reference  herein,  the  reflectivity  factors  computed  from  Eqs.  (2) 
and  (3)  will  be  designated  the  "reflectivity  factors  of  radar  measurement". 

The  radar  reflectivity  factors  may  also  be  computed  from  measured  or  as¬ 
sumed  knowledge  of  the  size  distribution  properties  of  the  hydrometeors.  If  the 
knowledge  stems  from  measurements,  the  reflectivity  factor  for  water  hvdro- 
meteors  is  given  by 


i=D 


max 


N  (D;)  Dj6  . 


(4) 


3 

where  N  (Dj)  is  the  number  of  the  drops  per  unit  volume  (per  m  )  of  the  classified, 

mid-diameter  size,  D.  (in  mm),  and  where  summation  is  made  from  the  minimum 

observed  diameter,  D  .  ,  to  the  maximum,  D  .  Analogously,  the  reflectivity 
min  m  sx 

factor  for  ice  hydrometeors  is  given  by 


i=D* 


min 


(5' 


where  De.  is  the  "equivalent  melted  diameter"  of  the  snow  or  ice  crystal  particlr 
of  the  ith  equivalent-diameter  class.  Assumptions  must  be  made  in  the  use  of 
Eq.  (5)  concerning  the  "effective  density"  of  the  snow  aggregates  or  ice  crystals. 
These  are  necessary  to  establish  the  relations  between  the  physical  sizes  of  the 
particles  and  their  equivalent  melted  diameters. 

The  subscript  "M"  applied  to  Zw  and  Z(  in  Eqs.  (4)  and  (5)  signifies  that 
these  are  particular  reflectivity  factors  computed  from  direct  measurements  of 
size  distribution. 

Values  of  the  radar  reflectivity  factors  may  additionally  be  obtained  by  use  of 

theoretical  distribution  functions  which  approximate  the  actual  distributions:  see, 

5  2  4 

for  example,  Laws  and  Parsons,'  Wexler,  Marshall  and  Palmer,  Boucher  and 
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Bartnoff,28  Gunn  and  Marshall, 14  .loss  et  al.29  and  Sekhor.  and  Srivastava.23  In 
such  case. 


D 

f  r 


max 


t  =  Jn 


N  [f(D)1  D°  dD 


min 


and 


,,  max 

N 

JDe 


(f(De)]  De°  dDe 


e 

min 


(6) 


(7) 


where  N  [f(Dl]  ar  ,i  De>]  are  particular  assumed  distribution  functions  and 
where  the  subscript  "T"  on  Z^y  and  Zj  indicate  the  theoretical  nature  of  the  re¬ 
flectivity  factors  thus  computed. 

To  illustrate  this  method  of  computation,  reference  is  made  to  the  work  of 
4 

Marshall  and  Palmer  concerning  rain.  They  assumed  a  distribution  function  of 
exponential  type, 


N  [f(D>]  =  Nd  =  a  e"AD  , 


(8) 


where  ND  is  the  number  of  drops  of  diameter  D  per  unit  volume  per  diameter 
bandwidth.  In  addition,  t.iey  assumed  that  the  factor  A  in  the  exponent  was  func¬ 
tionally  dependent  on  the  precipitation  rate,  P,  in  the  manner, 


A 


41  P 


-0.21 


-1 

cm 


(9) 


where  the  values  of  the  constant  and  exponent  were  deduced  from  observed  drop 
size  distributi  ins  of  rain. 

By  the  substitution  of  Eq.  (9)  into  (8)  into  (6),  with  integration  between  the 
diameter  limits  zero  to  infinity,'11 


28.  Boucher,  R..I.,  and  Bartnoff,  S.  (1955)  A  Comparison  of  Theoretically  Derivec 

and  Observed  Drip-Size  Distributions  in  Clouds  and  Hain.  Tufts  ['niversitv. 
Dept,  of  Physics,  Sci.  Upt  No.  4  under  Contract  AF  19(6  041-550,  30  pp. 

29.  .loss,.!.,  Schram,  Karin,  Thams,  .1.  C.  ,  and  Waldvogel,  A.  (1970)  On  the 

Quantitative  Determination  of  Precipitation  bv  Radar.  Wissenschaftliche 
Mitteilung  Nr.  63  I’orschungsstelle  der  "Fidgenossischen  Kommission  zum 
Studium  der  Ilagelbildung  und  der  Hagelabwehr"  am  Osservatorio  Ticinese 
della  Centrale  Meteorologica  Svizzera,  Locarno-Monti. 

*Drops  of  infinite  size  do  not  physically  exist,  of  course,  but.  the  equations  apply 
equally  well,  for  practical  purposes,  to  any  hydrometeor  population  in  which 
the  largest  drops  have  diameters  approximately  >  15/A  (see  Section  7). 
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Zw,r  =  3700  a  P1'47  mm6  m‘3  .  (10) 

4 

Marshall  and  Palmer  made  the  further  assumption  that  a  =  0.  08 ,  again 
supported  by  data  which  yielded  the  equation 

ZVVt  -  296  P1'47  mm6  m"3  ,  (ill 

which  was  referenced  earlier  in  the  present  report. 

The  Marshall  and  Palmer  method  indicated  above,  involving  the  use  of  a  dis¬ 
tribution  function  of  exponential  type,  has  become  a  more  or  less  standard  method 
for  determining  the  vs  P  relationships  for  rain,  also  the  '/.j^  vs  P  rela¬ 

tionships  for  snow  and  ice  crystals  (see  Gunn  and  Marshall,  Sekhon  and 

23  1 5 

Srivastava,  and  Carlson  and  Marshal!  '  ). 

Various  efforts  have  been  made  to  truncate  the  distribution  function  at  diam- 

23  30 

eter  limits  other  than  zero  and  infinity,  as  by  Sekhon  and  Srivastava,  *  but 
this  rapidly  leads  to  mathematical  complexity  and  requirements  for  additional 
assumptions. 

It  can  be  seen  that  the  distribution-function  method  requires  observational 
data  to  provide  realistic  estimates  of  the  constant  a,  and  of  the  coefficient  and  ex¬ 
ponent  of  the  Kq.  (9)  expression  (the  general  expression)  for  A  . 

3.2  Tlir  Hydromrleor  Parameters 

We  turn  now  to  a  more  detailed  consideration  of  two  of  the  hydromet’or  pa¬ 
rameters  of  SAMS  interest,  the  precipitation  rate,  P,  and  the  liquid-water-content, 
M  . 

The  precipitation  rate  is  a  volume  flux  parameter  that  tells  the  volume  of 

liquid  water  (in  the  form  of  drops  or  ice  particles'  that  falls  across  a  horizontal 

surface  per  unit  surface  area,  per  unit  time.  Thus  P  has  the  dimensions  LT  , 

and  it  is  conventional  to  employ  the  particular  units  mm  hr  '  . 

The  precipitation  liquid-water-content  is  a  density  parameter  that  specifies 

the  mass  amount  of  liquid  water  (of  precipitable  size)  that  exists  within  a  unit  at- 

-  3 

mospheric  volume.  "  The  dimensions  of  M  are  Mass  X  I,  ,  and  the  particular 
units  gm  m  3  are  .  ustomarily  utilized. 


30.  Sekhon,  R.S.  ,  and  Srivastava,  It.  C.  (1971)  Doppler  radar  observations  of 
drop  size  distributions  in  a  thunderstorm,  .1.  Atmos.  Sci.  28:983-994. 

-  n  ii 

A  distinction  is  made  between  the  liquid-water-content  of  precipitation-size  hydro¬ 
meteors,  generally  larger  than  about  80  microns  diameter  (or  equivalent 
melted  diameter),  which  can  be  detected  by  the  SAMS  radars,  and  the  liquid- 
water-content  of  cloud-size  hydrometeors  of  smaller  diameter,  which  carnot 
be  detected  by  the  radars  and  must  be  measured  by  aircraft  instruments. 
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The  precipitation  rate  and  liquid-water  content  may  be  related,  if  the  fall 
velocity  of  the  hydrometeors  can  be  specified  as  a  function  of  their  size.  Such 
specification  may  be  made  on  a  size-class  by  size-class  ba.-.is,  if  actual  fall  veloc¬ 
ity  data  are  utilized,  or  it  may  be  accomplished  on  the  basis  of  an  equation  rela¬ 
tion  that  describes  the  approximate  variation  of  fall  velocity  with  diameter  (equiv¬ 
alent  diameter  in  the  case  of  snow  or  ice  crystals).  One  commonly  used  relation 

3 

that  pertains  to  particular,  restricted  size  ranges,  seen  in  the  works  of  Spilhaus, 
32  18 

Gunn  and  Kinzer.'  and  Langleben,  for  example,  is 

V  =  a  Db  .  (12) 

for  water  hvdrometeors,  or 

V  =  a  H  b  .  (13) 

e 

for  ice  hydrometeors,  where  V  is  the  fall  velocity,  usually  in  cm  sec  *,  and  D 
and  De  are  the  diameter  and  equivalent-melted-diameter,  respectively. 

At  the  ground  surface,  it  is  usually  the  precipitation  rate  that  is  measured 
directly  (by  rain  gauges,  disdrometers,  nr  other  means)  and  the  liquid-water- 
content  is  derived  therefrom.  Above  the  ground  surface,  on  the  other  hand,  in 
storm  regions  where  aircraft  measurements  or  indirect  probing  methods  are  re¬ 
quired.  it  is  usually  the  liquid-water-content  that  is  measured  directly  and  the 
precipitation  rate  is  derived. 

3.  2.  1  MKASFKKMKNTS  AND  DERIVATIONS 
AT  Till]  GROPNI)  I.KVKI, 

For  example,  for  rain  at  the  ground  level  when  measurement  knowledge  exists 
concerning  the  size-distribution  flux  of  the  raindrops,  [as  acquired  bv  disdrom¬ 
eters,  for  instance  (see  Plank1  for  description)]  the  precipitation  rate  is  given  by 

i=  I) 

fiw  max  . 

P  =  ^  Nf(D.)  IF  mm  hr  ,  (14) 

i-l)  . 
min 

where  Nf(l>.)  is  the  number  of  drops  of  the  classified  diameter  si/e  I).  ,  in  mm, 

I  i  p  i 

that  fall  across  a  horizontal  sampling  area  A  ,  in  cm  ,  in  a  time  interval  At  in 
seconds. 


31.  Spilhaus,  A.  F.  (1048)  Raindrop  size,  shape,  and  falling  speed,  .1.  Meteorol. 
.5:108-1  10. 

Gunn,  R.  ,  and  Kinzer,  G.  I).  <10401  The  terminal  velocity  of  fall  for  water 
droplets  in  stagnant  air,  .1,  Meteorol.  6:243  (5h5,  504,  500-7). 
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The  liquid  water  content  of  this  rain  may  be  determined  from  information 
about  the  fall  velocity  of  the  drops  (from  actual  data  or  equation  relations,  as  men¬ 
tioned  earlier!. 

If  we  assume  that  Kq.(l2)  describes  the  fall  velocity  of  the  drops,  the  liquid 
water  content  is  determined  as  follows: 

The  wa'.or  mass  content  of  any  given  size-class  of  the  drops  is 

it  y  io‘3  p  n.(d.)  o3 

Massii  =  - 7. - - - —  S111  »  H5) 

i  6 

_3 

where  p  is  the  density  of  water  in  gm  cm  ,  and  1).  is  measured  in  mm. 

The  atmospheric  volume  (average  volume)  that  contained  these  drops  just 
prior  to  their  gravitational  fall  onto  the  collector  is  approximately, 


VolDi  =  VDj  As  At  •  06) 

where  V[).  is  the  fall  velocity  of  the  drops.  This  fall  velocity,  from  Eq,  (12),  is 
V.j  =  aIJ  b  .  (17) 

l  1 


The  liquid  water  content  of  the  drops  of  this  size  class  is,  by  definition,  the 
water  mass  per  unit  atmospheric  volume.  Hence,  from  Eqs.  (15)  and  (16),  also 
replacing  V|j.  in  Kq.  (16)  by  the  Eq.  ( 17)  expression, 
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where  Mi,  is  given  in  the  conventional  units  of  gm  m  (for  A  specified  in  cm  , 

i  s 

with  At  in  seconds). 

The  total  liquid  water  content  of  the  distribution  of  all  drop  diameters  is  then 
given  simply  by 
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The  precipitation  rate  at  the  surface  level  may,  of  course,  also  be  measured 
without  detailed  knowledge  of  the  size  distribution  of  the  raindrops.  Tipping  bucket 
or  weighing-tvpe  gauges  can  he  used  to  acquire  such  information. 
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For  P  values  obtained  directly  in  this  way,  it  has  been  customary,  in  the 
SAMS  work  at  AFCRL,  to  assume  that  the  precipitation  liquid-water-content  is 
related  to  P  as, 

M  =  k  Pt  .  (20) 


where  k  and  e  have  differing  values  dependent  on  rainfall  type.  The  particular 

values  used  will  be  cited  later.  Equations  of  this  form,  relating  M  and  P,  have 

4  14 

been  described  by  Marshall  and  Palmer,  Gunn  and  Marshall,  and  others. 

Although  it  is  possible  to  measure  the  precipitation  rate  at  the  surface  level 
in  snow  using  rain  gauges  that  heat  and  melt  the  frozen  particles  as  they  fall  into 
the  collector,  such  instruments  were  not  used  at  Wallops  Island.  All  of  the  mis¬ 
sile  firings  thus  far  were  made  into  storms  having  rain  at  the  surface  level. 

3.  2.  2  MEASUREMENTS  AND  DERIVATIONS  ALOFT 

At  storm  levels  above  the  surface  level,  aircraft,  or  other,  measurements  of 
hydrometeors  usually  consist  of  volumetric  samples  that  provide  information  about 
(1)  the  types,  sizes,  and  numbers  of  the  drops  or  particles,  or  (2)  the  bulk  or  in¬ 
tegrated  properties  of  the  hydrometeors,  such  as  their  liquid  water  content. 

If  size  distribution  measurements  are  obtained  from  aircraft  by  the  use  of 
spectrometer,  formvar  replicator,  or  foil-sampling  devices,  for  example,  the 
liquid  water  content  of  the  distribution  is  given  by. 
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(21) 


where  N  (I).)  is  the  number  of  the  drops  of  the  classified  diameter  size  D.  (or 
8  1  r  1 

Dp.  ,  in  the  case  of  ice  hydrometeors)  in  the  sampling  volume,  V  ,  which  latter  is 
“i  3  s 

specified  in  m  . 

The  precipitation  rate  may  be  computed  from  these  aircraft,  size-distribution 
measurements  by  utilizing  the  fall-velocity  data  or  relations  cited  earlier  and  by 
making  the  assumption  that  no  updraft  or  downdraft  motions  exist  in  the  storm  at 
the  altitude  of  interest.  If,  for  example,  we  use  Eq.(12)  to  specify  the  fall  veloc¬ 
ities  of  the  hydrometeors,  the  precipitation  rate,  following  the  same  reasoning 
steps  described  in  Section  3.  2.  1  ,  except  in  reverse,  is 
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When  the  liquid-water-content  is  measured  directly  by  the  aircraft,  as  by 
means  of  Australian*  or  total  liquid  water  content  meters,  or  certain  spectrometer 
devices,  the  corresponding  precipitation  rate  is  inferred  from  the  inverted  form  of 
Eq.  (20),  that  is 


(23) 


where,  as  previously  mentioned,  k  and  e  have  particular  values  dependent  on  the 
rain,  snow,  or  ice  crystal  typos  detected  by  the  aircraft.  The  assumption  of  a 
static  atmosphere,  devoid  of  updraft-downdraft  motions,  is  also  implicit  in  this 
equation,  as  with  Eq.  (22). 

3.3  Pr  viou*  Correlation*  and  Sampling  Volume  Contideration* 

Thus  far,  we  have  discussed  the  background  of  the  Z  vs  P  relationships, 

have  indicated  the  several  methods  of  computing  Z,  and  have  described  how  P 

and  M  are  measured  or  derived  at  the  ground  level  and  aloft.  Comment  will  now 

be  made  about  the  nature  of  the  correlation  and  regression-analyses  performed  in 

the  past  and  about  some  of  the  problems  involving  comparison  sampling  volumes. 

Much  of  the  previous  correlation  work  which  yielded  particular  regression  or 

best-fit  equations  of  Z  vs  P,  as  summarized  in  Section  2,  was  performed  using 

Z  and  P  values,  both  of  which  were  computed  from  measurements  (or  presumed 

knowledge)  of  the  size  distribution  properties  of  the  hydrometeors.  In  other  words, 

the  Z  and  P  values  used  to  obtain  the  regression  or  best-fit  equations  [of  the 

Eq.  (1)  type  herein]  were  not  independent  values,  rather  they  were  interrelated 

through  the  commonality  of  the  size-distribution  data  used  for  computations.  Ad- 

33 

ditionally,  as  discussed  by  Mason,  many  of  these  equations  were  obtained  from 

data  acquired  from  sampling  volumes  of  relatively  small  size,  of  the  order  of  0.1 

to  1.0  m* .  This,  Mason  notes,  probably  contributes  to  the  considerable  differ- 

34 

enres  among  the  equations  reported  in  the  literature.  Mueller  and  Sims  have 
also  discussed  this  problem. 


*The  Australian  liquid-water-content  meter  is  an  instrument  in  which  special 

electrolytic  paper,  on  a  roll,  is  exposed  to  the  uirstream,  behind  a  slit.  The 
impinging  water  droplets  moisten  the  paper  and  change  its  conductivity,  pro¬ 
viding  a  measure  of  the  liquid-water-content. 

33.  Mason,  B.  .1.  (1071)  The  Physics  of  Clouds,  Second  Edition,  Clarendon  Press, 

Oxford,  England. 

34.  Mueller,  E.A.,  and  Sims.  A.  L.  (1066)  The  influence  of  sampling  volume  on 

raindrop  size  spectra,  Proc.  12th  Weather  Radar  Conf. ,  Norman,  Okla¬ 
homa,  p.  135(470,  6001. 


21 


The  difficulty  with  small  sampling  volumes  is  that  the  largest  hydrometeors, 

3 

which  exist  in  relatively  small  number  concentration  in  the  atmosphere  (1  Der  m 
3 

to  1  per  100  m  ,  or  so)  are  either  not  represented  or  are  inadequately  repre¬ 
sented  in  the  sample.  These  largest  hydrometeors  do  not  contribute  significantly 

to  the  liquid-water-content  or  precipitation  rate,  which  are  essentially  functions 
3 

of  D  ,  D  being  diameter,  but  they  do  contribute  importantly  to  the  reflectivity 

6 

factor  which  is  dependent  on  D  .  In  fact,  a  few  large  hydrometeors  may  con¬ 
tribute  more  to  the  reflectivity  factor  than  does  the  totality  of  all  the  other  smaller 
hydrometeors.  Thus,  one  must  be  exceedingly  cautious,  in  the  use  of  empirical 
Z  vs  P  relationships,  to  insure  that  these  were  obtained  from  data  samples  of 
sufficient  size  to  be  representative  of  the  pulse  volume  (or  integration  volume)  of 
the  radar. 

It  would  seem  that  the  Z  vs  P  relationships  of  greatest  validity  would  hr 
those  obtained  from  correlation  analyses  in  which  radar- measured  values  of  Z 
|fora  calibrated  radar,  see  Eqs.  (2)  and  (3)]  were  correlated  with  independent 
values  of  P,  acquired  from  surface  or  aircraft  measurements.  Such  studies,  in- 

3 

volving  surface  rainfall  measurements,  have  been  reported  by  Marshall  et  al, 

35  36  37  38 

Hooper  and  Kippax,  Hood,  Austin  and  Williams,'  and  Doherty.'  These 

direct  correlation  investigations,  however,  have  been  relatively  few  in  number, 

primarily  because  of  the  difficulties  and  uncertainties  of  radar  calibration  and 

39 

proper  pulse-volume-averaging,  as  noted  by  Austin.  Questions  exist  about  the 
validity  of  various  of  the  previous  results. 

I.  THE  AKCRL  MEASUREMENT  AND  ANALYTICAL  PROGRAM 

The  objectives  of  the  AFCHL  measurement-analytical  program  for  the  SAMS, 
Wallops  Island  storms  may  now  be  described. 


35.  Hooper,  .1.  E.  N.  ,  and  Kipp  .A. A.  (1950)  Radar  echoes  from  meteorological 

precipitation,  Proc,  I,  E.  ,  97 ( Pt.  1):B9. 

36.  Hood,  A.l).  (1950)  Quantitative  measurements  at  3  and  10  centimetres  of 

radar  intensities  from  precipitation,  Nat,  Res.  Counc.  ,  Rpt  No.  2155, 
Ottawa. 

37.  Austin,  P.  M.  ,  and  Williams,  E.I..  (1951)  Comparison  of  radar  signal  inten¬ 

sity  with  precipitation  rate,  M.  I.  T.  Weather  Radar  Research,  Tech.  Rpt 
No'.  14. 

38.  Doherty,  I..  II.  (1963)  The  scattering  coefficient  of  rain  from  forward  scatter 

measurements,  Proc.  of  the  Tenth  Weather  Radar  Conf.  ,  A.M.S. ,  Boston, 
p.  171. 

39.  Austin,  I’.  M.  ( 1  *^64 )  Radar  measurements  of  precipitation  rate,  Proc.  of  the 

Eleventh  Raoar  Conf. 


22 


The  general  objectives  were  to  obtain  observational  and  measurement  infor¬ 
mation  for  each  of  the  storms  that  would  (1)  define  the  hydrometeor  types  and 
regions  that  were  present  in  the  storms  at  the  different  levels,  (2)  establish  the 
particular  relationships  among  Z,  P,  and  M  that  pertained  to  the  observed  hydro¬ 
meteor  types,  (3)  utilize  these  relationships  to  compute  the  M  and  P  values  cor¬ 
responding  to  the  radar- measured  values  of  Z  for  the  missile  trajectories,  and 
(4)  also  compute,  or  otherwise  determine  the  size-distribution  properties  of  the 
hydrometeors  along  the  trajectories. 

4.1  The  llydronieteor  Regions 

There  are  four  identifiable  hydrometeor  regions  that  exist  in  the  winter 
storms  at  Wallops  Island  which  may  be  defined  in  terms  of  the  type  and  size- 
distribution  properties  of  the  hydrometeors.  These  regions,  from  tep  to  bottom 
of  the  storm,  in  the  order  of  the  progressive  development  of  the  hydrometeors, 
are  (1)  the  ice-crystal  region,  (2)  the  small-snow  region,  (3)  the  large-snow 
region,  and  (4)  the  rain  region.  The  hydrometeors  in  each  of  these  regions  may, 
in  turn,  be  sub-classified  on  the  basis  of  their  particular  characteristics.  Thus, 
sub-classification  in  the  ice-crystal  region  is  made  on  the  basis  of  the  specific 
type  of  ice-crystals  that  are  present  as,  for  example,  columns,  plates,  dendrites, 
or  needles.  In  the  snow  regions,  it  is  dependent  on  the  predominant  types  of  ice- 
crystals  that  are  contained  in  and  comprise  the  snow  aggregates.  In  the  rain  re¬ 
gion,  it  is  made  on  the  basis  of  the  nature  of  the  size-distribution  of  the  raindrops. 

The  hydrometeor  regions  and  sub-classes  used  in  the  SAMS  program  are 
identified  in  Table  2. 

Not  all  of  the  hydrometeor  regions  are  necessarily  present  in  any  given  storm. 
If  the  storm  has  a  convective  nature,  for  example,  and  depending  on  how  the  SAMS 
missile  traversed  the  storm,  there  may  not  be  an  ice  crystal  region  at  the  top  of 
a  particular  convective-element  intercepted  by  the  missile,  or  perhaps  neither  the 
ice  crystal  region  nor  small-snow  region  may  exist  in  such  a  situation. 


4.2  The  Kqualion*  l  urd 

Particular  empirical  equations  of  the  types  described  in  Sections  2  and  3 
were  used  for  each  of  the  above  hydrometeor  regions  to  determine  the  values  of 
P  and  M  that  corresponded  to  the  radar-measured  values  of  Z.  Thus,  P  was 
computed  from  the  equation 


P 


K 


P 


(24) 


which  is  an  inversion  of  the  traditional  form  of  relationship  of  Kq.  (1). 
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•A  type  used  only  in  the  1971-72  season. 


The  Kp  and  Ep  values  of  this  equation  were  specified  for  each  of  the  hydro¬ 
meteor  regions  and  each  hvdrometeor  type  (where  information  existed).  The  par¬ 
ticular  values  used  for  the  data  of  the  1971-72  and  1972-73  SAMS  seasons  are 
listed  in  Table  2.  The  values  for  the  rain  region  were  obtained  from  disdrometer 
data  and  from  the  work  of  .loss  et  al.^®  The  values  for  the  snow  and  ice  crystal 
regions  were  obtained  by  Cunningham  from  analyses  that  will  be  described  in  a 
subsequent  report. 

The  Z  values  employed  in  the  above  equations,  except  for  the  surface  and 
uppermost  layers  of  the  storm  (these  exceptions  will  be  discussed),  were  the 
radar- measured  values  of  Z w,  as  specified  by  Eq.  (2),  for  water  hydrometeors, 
and  Zj ,  as  specified  by  Eq.  (3),  for  ice  hydrometeors. 

It  was  assumed,  in  all  of  the  storm  ana.vses,  that  the  precipitation  liquid- 
water-content  was  re'nted  to  the  precipitation  rate  in  the  manner  described  by 
Eq.  (20)  which,  rewritten,  is 

M  =  kPe  ,  (20) 

whore  again  the  particular  values  of  k  and  t  were  specified  for  each  hydrometeor 
region  and  type,  as  shown  in  Table  2. 

It  follows,  from  Eqs.(24'  and  (201,  that  the  precipitation  liquid  water  content 
is  related,  or  implicitly  assumed  to  be  related,  to  the  radar  Z  values  as 

EM 

M  =  Z  IVI  ,  (251 

where 

K„  =  kK(  (261 

M  p 

and 

ems£V  (27) 

The  values  of  and  E^  for  E q .  (25)  were  established  for  each  hydrometeor 

region  and  type  either  (1)  from  knowledge  of  K  ,  E  ,  k  and  e,  using  Eqs.  (26) 

^  P 

and  (27),  or  (2)  from  regression  analyses  in  which  aircraft-determined  values  of 
\1  and  Z,  obtained  from  size  distribution  measurements,  were  inter-correlated. 

In  either  event,  the  consistency  of  the  equation  set,  composed  of  Eqs.  (24),  (20), 


40,  .loss,.!.,  Thams,  I.C. ,  and  Waldvogel,  A.  (1968)  The  variation  of  raindrop 

size  distributions  at  Locarno,  Proc.  Into rnatl.  Conf,  on  Cloud  Physics, 
Toronto,  Amer.  Meteorol.  Soc. ,  Boston,  p.  369. 
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and  (25),  was  maintained  through  use  of  the  Eqs.(26)  and  (27).  The  details  of 
these  analyses  will  be  described  by  Cunningham,  as  mentioned  earlier.  The  par¬ 
ticular  K^j  and  values  determined  from  his  analyses  are  shown  in  Table  2. 

It  will  be  noted,  with  reference  to  Table  2,  that  the  coefficients  and  exponents 
are  presented  for  only  one  of  the  ice  crystal  types  listed  and  for  only  one  of  the 
small-snow  types.  This  reflects  the  fact  that  AFCRL  has  not  yet  established  the 
values  for  the  other  listed  types.  One  of  the  prime  objectives  of  the  continuing 
SAMS  program  at  AFCRL  is  the  acquisition  of  appropriate,  simultaneous  radar 
and  aircraft  information  that  will  permit  the  specification  of  the  equation  values 
for  these  other  hydrometeor  types  (and  that  will  also  permit  the  accuracy  verifica¬ 
tion  of  all  equations  for  all  hydrometeor  types). 

4.3  The  Equations  in  Inverse  Form 

Most  of  the  computations  for  the  SAMS  missile  trajectories  were  performed 
using  Eqs.(24),  (20),  and  (25).  However,  occasional  need  also  existed  for  com¬ 
putations  employing  the  inverse  form  of  the  equations,  in  which  the  "independent" 
and  "dependent"  variables  are  reversed.  The  constants  and  expoi  n;s  of  the 
equations  in  this  form  are  presented  in  Table  3.  These  are  the  more  or  less  tra¬ 
ditional  forms  of  the  equations;  hence  their  presentation  facilitates  comparisons 
with  the  equations  of  other  investigators,  as  reported  in  the  literature. 

4.4  Reflectivity  Factor  Values  for  the  Missile  Trajectories 

The  radar  values  of  the  integration  signal,  I ,  were  determined  for  the  major 
portions  of  the  SAMS  missile  trajectories  by  methods  described  in  the  previous 
report  by  Plank.'  (Subsequently,  this  report  will  be  referred  to  as  "R  No.  1".) 
Special  problems  existed,  however,  in  determining  or  estimating  the  I  values 
that  existed  within  two  particular  altitude  layers  within  the  storms,  namely  the 
surface  layer  and  the  highest  layer. 

Ground  clutter  return  from  the  launch  site  location  masked  the  hvdrometeor 
echoes  along  the  lowest  portion  of  the  missile  trajectories  (to  about  0.4  km  alti¬ 
tude),  and  it  was  necessary  to  assume  that  the  I  signals  received  from  just  off¬ 
shore  of  the  launch  site  (about  1.3  km  offshore)  were  representative  of  those  of 
the  surface  layer  trajectories  at  the  same  altitudes. 

The  backscatter  return  from  the  hvdrometeors  along  the  trajectory  portions 
near  and  at  the  storm  top  at  the  range  extremity  of  the  missile  paths  through  the 

-For  •  xample,  the  radar  values  of  Z  could  not  be  obtained  for  the  surface  level 
pow;t  of  the  missile  trajectory,  because  of  ground  clutter.  Thus,  for  purposes 
of  overall,  computational  consistency,  the  "effective  values  of  Z"  had  to  be 
determined  from  'lie  surfa"e  measurements  of  P,  which  involved  use  of  the 
"reverse-form  equations"  cited  above. 
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Table  3.  Inverse  Form  of  the  Equations  of  Table  2 


wigs 


»  rv  F'-r  •w,vrr,r«"  - 


27 


*A  type  used  only  in  the  1971-7  2  SAMS  season. 


storms  were  sometimes,  in  certain  storms,  below  the  minimum  detectable  by  the 
radar,  in  such  cases,  the  I  values  for  the  upper  trajectory  portions,  that  is,  for 
the  uppermost  layer  of  the  storms,  had  to  be  estimated.  For  storms  of  homo¬ 
geneous  structure,  the  assumption  was  made  that  the  I  values  detected  bv  the 
radar  at  vertical  incidence  would  be  representative  of  those  of  the  missile  tra¬ 
jectory  in  its  upper  portion.  *  For  storms  of  convective  structure,  on  the  other 
hand,  such  an  assumption  could  not  be  made,  because  of  the  spatial  variability  of 
the  radar  i  values  in  the  top  part  of  the  storm.  The  trajectory  values,  in  this  lat¬ 
ter  case,  could  merely  be  considered  to  be  something  smaller  than  the  minimum 
detectable  values  at  the  particular  range. 

Aircraft  observations  and  measurements  were  made  in  the  Wallops  storms 
during  the  periods  immediately  following  the  missile  firings.  These  observations 
and  measurements  provided  information  ( 1  >  about  the  altitude  boundaries  and  ver¬ 
tical  extent  of  the  melting  zone,  the  hydrometeor  regions  and  the  transition  zones 
within  the  storms,  and  (2)  about  the  types  and  size-distribution  properties  of  the 
hydrometeors  in  each  region  and  zone. 

Information  was  additionally  available  concerning  the  precipitation  rate  and 
size-distribution  of  the  raindrops  at  the  surface  level  (acquired  from  rain  gauges 
and  disdrometers)  and  concerning  the  vertical  temperature  and  humidity  structure 
of  the  storm  (acquired  from  aircraft  and  special  radiosonde  measurements). 

The  aircraft  information  about  the  altitude  boundaries  of  the  melting  zone  per¬ 
mitted  the  differentiation  of  the  regions  of  water  and  ice  hydrometeors.  This,  in 
turn,  permitted  the  reflectivity  factors,  Z^y  and  Zj ,  to  be  computed  from  the 
trajectory  values  of  the  radar  integration  signal,  I .  The  procedures  and  equa¬ 
tions  used  have  been  described  in  R  No.  1.  It  was  also  noted  in  that  report  that 
both  of  the  reflectivity  factors,  Z^, ,  for  water,  and  Zj ,  for  ice,  were  computed 
across  the  melting  zone  itself. 

4.5  Computation*  of  Precipitation  Rate  and  Liquid  Water  Content 

The  trajectory  values  of  P  and  M  were  computed  from  the  Z^y  and  Zj 
values  using  the  methods  and  assumptions  described  in  the  following  sections. 

The  description  begins  at  the  storm  top  level  and  continues,  section  by  section 
within  the  text,  for  each  of  the  hydrometeor  regions  and  zones  that  is  normally 
observed  in  a  Wallops  Island  storm  of  homogeneous  (non-conveetive)  type.  Thus 
the  descriptive  order  follows  that  of  the  natural  developmental  processes  that 


*The  radar  at  vertical  incidence  can  detect  the  I  signals  from  the  storm  top 

whereas  the  radar  at  the  elevation  angle  of  the  missile  trajectory  cannot.  The 
reason  is  that  the  radar  range  t<>  the  storm  top  >s  smaller  at  vertical  incidence. 
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occur  in  the  storm  as  the  hydrometeors  are  generated  and  fall  gravitationally 
through  the  depth  of  the  storm  to  the  surface  level. 

Comments  about  the  observed  nature  of  the  regions  and  zones,  and/or  about 
particular  computational  problems,  are  made  at  the  beginning  of  each  section. 

The  computation  procedures  used  to  obtain  P  and  M  art  explained  thereafter. 

4.  5.  1  STORM  TOP 

The  clouds  at  the  upper  boundaries  of  the  storms,  at  the  exit  points  of  the 
SAMS  missiles,  were  not  reliably  detected  by  the  radar,  at  least  not  in  the  1970 
to  1973  seasons.  Thus,  the  storm-top  altitudes  at  the  .nissile  exit  points  had  to 
be  established  either  from  aircraft  observations,  when  aircraft,  usually  the  AFCRL 
C-130,  were  present  that  could  fly  at  these  altitudes  (of  approximately  26,000  to 
32,  000  ft),  or  they  had  to  be  estimated  from  vertical-incidence  radar  data  and/or 
from  the  humidity  profiles  or  radiosonde  ascents.  Estimation  accuracy,  in  the 
latter  case,  was  probably  about  ±  1000  ft  or  so  for  storms  with  uniform  top  struc¬ 
ture,  and  was  possibly  no  better  than  ±2000  ft  to  ±5000  ft  for  storms  having  con¬ 
vective  top  structure,  dependent  on  the  spatial  organization  and  altitude  variability 
of  the  convective  elements. 

The  radar  ZW  value,  also  the  P  and  M  values,  were  assumed  to  have  zero 
values  at  the  storm  top  level.  This  may  be  patently  obvious,  but  is  mentioned  for 
completeness. 

4.  5.  2  ICE  CRYSTAL  REGION 

The  uppermost  portion  of  most  of  the  Wallops  storms  (again  excepting  those 

of  convective  type)  is  comprised  of  populations  of  individual  ice  crystals  which 

exist  at  altitudes  above  about  25,000  ft  to  28,  000  ft,  where  the  temperatures  are 

generally  less  than  -25°C.  These  ice  crystals  are  generated  and  developed  in  a 

supersaturated  environment  (with  respect  to  ice)  which  is  created  by  the  uplift  mo- 

41  42  43 

tions  within  the  storm  in  general  (see  Weickmann,  Magono,  Borovikov, 

44  45  46 

Atlas,  Hobbs,  and  Ono).  The  ice  crystals  have  various  different  shapes, 

41.  Weickmann,  H.  (1947  Die  Eisphase  in  der  Atmosphare.  Library  Trans.  273, 

Royal  Aircraft  Establishment,  Farnborough,  96  pp. 

42.  Magono,  C.  (1953)  On  the  growth  of  snowflakes  and  graupel. ,  Scient.  Rpt 

Yokohama  Univ.  ,  Ser.  1,  No.  2,  p.  18(239,  24 1,  242). 

43.  Borovikov,  A.  M.  (1953)  Some  results  on  an  investigation  of  the  structure  of 

crystal  clouds  (in  Russian)  Trudy  Tsentral,  Acrolog.  Obs.  No.  12(247). 

44.  Hosier,  C.  L.  ,  Jensen,  D.  C. ,  and  Goldshlak,  L.  (1957)  On  the  aggregation 

of  ice  crystals  to  form  snow,  J.  Meteorol.  14:415(250). 

45.  Hobbs,  P.  V.  (1969)  Ice  multiplication  in  clouds,  J.  Atmos.  Sci.  26:315-318. 

46.  Ono,  A.  (1970)  Growth  mode  of  ice  crystals  in  natural  clouds,  J.Appl.  Sci. 

27(4):649-658. 
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dependent  on  the  temperature-humidity  conditions  of  their  formation.  Numerous 

experimental  investigations  have  been  conducted  to  ascertain  the  specific  crystal 

types  that  occur  under  different  meteorological  conditions;  for  example,  see 
47  4ft  4Q  ftO  ft  1  ft? 

Ludlam,  Nakaya,  Gold  and  Power,  *'  u  Hosier,  Magono  and  Lee,  and 

Kikuchi.^ 

The  generally-observed,  most-common  types  of  ice  crystals  are  identified  in 
Table  2. 

The  P  and  M  values  for  those  portions  of  the  missile  trajectories  that  passed 
through  the  ice  crystal  region  were  computed  from  the  radar  Zj  values  using  the 
equations  shown  in  Table  2.  Thus,  P  waa  first  computed  from  Zj ,  employing  the 
equation  of  the  first  column,  and  M  was  then  computed  from  P,  utilizing  the 
equation  of  the  second  column.  The  reason  for  computing  P  first  is  related  to 
the  transition-zone,  interpolation  problem  discussed  in  the  following  section. 
Actually,  however,  the  computational  order  is  unimportant,  since  the  equations  of 
Table  2  are  self  consistent  for  each  of  the  hydrometeor  types. 

Only  one  set  of  equations  for  ice  crystals  is  presented  in  Table  2,  which  per¬ 
tains  specifically  to  ice  crystals  consisting  of  columns,  bullets,  and  rosettes. 
These  equations,  as  mentioned  earlier,  are  the  only  ones  presently  available. 
Consequently,  they  had  to  be  used  commonly  for  all  of  the  trajectory  computations 
of  P  and  M  in  the  ice  crystal  region,  irrespective  of  the  actual  crystal  types  ob- 
serveu  ir  the  Wallops  storms. 

4.  5.  3  ICE-CRYSTAL  TO  SMALL-SNOW  TRANSITION  ZONE 

As  the  largest  ice  crystals  in  the  uppermost  part  of  the  storm  (those  which 
have  greater  fall  velocities  than  the  other  smaller,  more-numerous  crystals  of 
the  populations)  descend  through  the  supersaturated  environment  of  the  storm, 
they  grow,  presumably  both  by  diffusion  and  by  colliding  with  and  "cohesively 
capturing  various  of  the  other  smaller  ice-crystals  along  their  fall  trajectories. 


47.  Ludlam,  F.  H.  (1947)  The  forms  of  ice  clouds,  Quart.  J.  Roy.  Meteorol.  Soc. 

74:39-56.  . 

48.  Nakaya,  . I.  (1951)  The  formation  of  ice  crystals.  Compendium  of  Meteorol.  . 

Boston,  Amer.  Meteorol.  Soc.  207-220, 

49.  Gold,  L.  W. ,  and  Power,  B.  A.  (1952)  Correlation  of  snow-crystal  type  with 

estimated  temperature  of  formation,  J,  Meteorol.  9:447  (246-7). 

50.  Gild,  L.  W. ,  and  Power,  B.  A,  (1954)  Dependence  of  the  forms  of  natural 

snow  crystals  on  meteorological  conditions,  .1.  Meteorol.  1  1:35  (247). 

51.  Hosier,  C.  L.  (1954)  Factors  governing  the  temperature  of  ice-crystal  for¬ 

mation  in  clouds,  Proc,  Toronto  Met.  Conf.  ,  1953,  p.253,  R.  Met.  Soc.  , 
London  (157). 

52.  Magono,  C. ,  and  Lee,  C.  W.  (1966)  Meteorological  classification  of  natural 

snow  crystu’s,  .1.  Fac.  Sci. ,  Hokkaido  Univ.  2(Ser.  VII):32.-335. 
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The  precise  growth  mechanisms  are  not  clearly  known,  but  aircraft  measurements 
and  observations  indicate  that  the  largest  ice  crystals  do  progressively  grow  and 
develop  into  "agglomerates ",  or  "small  snow"  particles.  (Theoretical  and  experi¬ 
mental  information  concerning  aggregation  processes  have  been  reported  by 
19  5'^  34  44  55 

Magono,  "  '  Marshall  and  l.angleben,'  Hosier  et  al,  Nakaya  et  al,'  Hosier 

and  llallogren, ‘  ®  and  Podzi mek.'"’^! 

The  altitude  zone  in  the  storm  in  which  this  aggregation  process  is  most  pro¬ 
nounced  (as  observed  bv  aircraft  (see  footnote  belowll  is  herein  defined  to  be  "the 
ice-crystal  to  small-snow  transition  zone".  Seemingly,  this  is  a  physically-real 
zone,  which  separates  the  region  of  ice  crystals  aloft  from  the  small-snow  region 
below.  The  zone  is  bounded  at  the  top  at  some  upper  level,  where  small  single 
ice-crystals  predominate  within  the  populations,  and  only  an  occasional  aggregate 
is  observed  which  contributes  insignificantly  to  the  liquid-water-content  of  the 
populations.  The  zone  extends  downward  to  some  lower  level,  where  aggregates 
are  observed  to  be  present  in  numbers  such  that  they  are  the  prime  contributors 
to  tile  liquid-water-oontent  of  the  populations. 

The  precipitation  rate  (PI  and  liquid- water- content  (M)  values  for  the  missile 
trajectory  segments  crossing  this  "ice-crystal  to  small-snow  transition  zone" 
were  computed  as  follows: 

The  precipitation  rate  was  assumed  to  be  linearly  var.  able  through  the  zone, 
increasing  downward.  Thus, 


,Jd  = 


*VH+ 


(28) 
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54.  Marshall,  .1.  S.  ,  and  l.angleben,  M.  P.  (1954)  A  theory  of  snow-crystal  habit 

and  growth,  .1.  Meteorol.  11:104(254-6). 
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55.  Nakaya,  U. ,  llanajima,  M.  ,  and  Muguruma,  .1.  (1958)  Physical  investigations 

on  the  growth  of  snow  crystals,  .1.  Fac.  Sci. ,  Hokkaido  Univ. ,  Ser.  IF,  5,  87 
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56.  Hosier,  C.  I..  ,  and  Ilallgren,  It.  K.  (1961)  Ice  crystal  aggregation,  Nubila  4: 

13(249). 

57.  Podzimek,  .1.  (1968)  Aerodynamic  conditions  of  ice  crystal  aggregation, 

Proc.  Internatl.  Conf,  Cloud  Phys.  ,  Toronto,  295-299. 

The  presence  of  aggregates  in  a  population  of  ice  crystals  can  be  detected  from 
aircraft  (1)  by  seeing  (hem  occasionally  strike  the  "snow-stick"  instrument 
aboard  the  aircraft,  (2)  by  seeing  them  "pass  by"  a  black  background,  such 
as  an  engine  nacelle,  or  (3)  by  detecting  the  occasional  presence  of  "larger 
than  ordinary"  particles  using  the  "raindrop  spectrometer"  instrument  of  the 
Al  CHI.  C-  1 30 A  aircraft. 

These  statements  may  be  clarified  somewhat  by  pointing  out  that  in  almost  all 
hvdrometeor  populations  (composed  of  ice-crystals,  or  of  snow,  or  of  rain), 
the  smallest  particles  or  drops  are  the  most  numerous  ones;  see  the  size- 
distribution  discussion  of  Section  7.  'Inus,  the  zone-boundary  definitions  above 
are  "verbally  awkward"  because  they  refer  to  the  "third  moment  properties" 
of  the  size  spectra  of  the  hydrometeors.  If  the  definitions  are  bothersome,  it 
is  suggested  that  the  render  turn  directly  to  Section  7. 
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where  is  the  zone  precipitation  rate  at  the  distance  d  below  the  base  of  the 
ice  crystal  region,  where  Pj,  is  the  precipitation  rate  at  the  base  of  the  ice 
crystal  region  and 


(79) 


where  P^s-p  is  the  precipitation  rate  at  the  top  of  the  small-snow  region  (com¬ 
puted  as  indicated  in  the  following  section)  and  is  the  total  vertical  depth  of 
the  transition  zone.  The  subscript  "3"  on  a  and  H  indicates  that  this  transition 
zone  is  the  third  such  zone  above  the  surface  level,  as  will  be  explained. 

This  assumption  implies  that  the  volume  flux  of  water  mass  across  horizontal 
surfaces  within  the  transition  zone  is  increasing  linearly  with  the  fall  distance  of 
the  hvdrometeor  particles.  This  is  consistent  with  the  progressive,  downward 
development  of  ihe  precipitation  process  within  the  storm. 

The  liquid- water-content  values  ir  this  transition  zone  were  computed  from 
the  Pd  values  using  the  equation 


(30) 


where  \ld  is  the  liquid- water-content  at  the  distance  d  below  the  base  of  the  ice 
crystal  region  and  where  k(,  and  are  the  coefficient  and  exponent  of  Eq.(20). 
which  have  the  particular  values  pertaining  to  ice  crystals,  as  listed  in  Table  2, 
column  2.  The  factors  j3^  and  y^  are  given  by 


0 


3 


and 
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where  kSg  and  are  the  coefficient  and  exponent  of  Eq.  (20'  that  pertain  to 
small-snow,  as  also  specified  in  Table  2,  column  2. 

*The  assumption  that  P  is  linearly  variable  with  altitude  is  considered  to  be  more 
physically  realistic  than  the  alternate  possible  assumption,  which  is  that  M 
is  linearly  variable.  This  is  because  the  liquid-water-content  in  storms  may 
be  "held  in  storage"’,  at  any  given  level  by  two  mechanisms:  (1)  by  the  differ¬ 
ent  fall  velocities  of  the  constituents  of  hydrometeor  populations  of  different 
type  and  (2)  by  the  presence  of  updraft-downdraft  motions  within  the  storms. 
An  assumption  that  P  varies  linearly,  on  the  other  hand,  since  P  is  a  flux 
parameter,  onlv  involves  the  second  mechanism.  The  detailed  reasons  are 
left  unstated. 


Thus,  the  values  of  M^  in  the  transition  zone  were  computed  under  the  as¬ 
sumption  of  the  linear  altitude  variation  of  k  and  e  of  Eq.  (20).  * 

4.5.4  SMALL-SNOW  REGION 

The  small-snow  region,  as  the  name  indicates,  is  defined  primarily  on  the 
basis  of  the  size  of  the  observed  snow  particles.  This  region  extends  downward 
from  the  base  of  the  transition  zone,  described  in  the  previous  section,  to  some 
lower  level  where  aircraft  measurements  indicate  that  the  equivalent  melted 
diameter  of  the  largest  snow  aggregates  is  about  0.5  mm  (which  corresponds  to 
physical  sizes,  dependent  on  snow  type,  ranging  from  approximately  1  to  3  mm). 

Thus,  small-snow  and  the  small-snow  region  are  defined  in  terms  analogous 
to  those  used  to  define  "drizzle"  in  the  case  of  water  hydrometeors.  The  maxi¬ 
mum  equivalent- melted-diameter  specified  for  small-snow  is  the  same  as  the 

maximum  actual  diameter  for  drizzle  droplets,  as  drizzle  is  conventionally  de- 
5  8 

fined;  see  Huschke,  for  example. 

The  P  and  M  values  fo.-  the  trajectory  segments  that  passed  through  the 
small-snow  region  were  computed  from  the  equations  of  the  first  and  second 
columns  of  Table  2.  The  only  equations  available  for  small-snow  at  present  are 
those  of  Marshall  and  Gunn10  as  modified  by  Cunningham  (see  Tables  2  and  3). 
These  had  to  be  used  commonly  for  all  storm  computations  of  the  past  SAMS 
seasons. 

4.5.5  SMALI.-SNOW  TO  LARGE-SNOW  TRANSITION  ZONE 

The  distinction  between  the  small-snow  and  large-snow  regions  in  the  Wallops 
storms  was  made  on  rather  tentative  and  hypothetical  grounds.  Likewise,  the 
definition  of  the  transition  zone  between  these  regions  was  provisional  and  con¬ 
jectural. 

Observations  and  measurements  obtained  by  the  AFCRL  C-130A  aircraft  in 
the  Wallops  storm  of  2  February  1973  suggest  that  some  "rapid  growth  process" 
occurred  in  the  populations  of  falling  snowflakes  when  the  particles  had  attained 
maximum  sizes  of  about  1  to  3  mm  (about  0.5  mm  equivalent  melted  diameter). 
When  the  aircraft  flew  at  levels  below  the  altitude  where  these  sizes  occurred,  it 
was  observed  that  the  maximum-size  particles  seemingly  became  "much  larger" 


58.  Huschke,  R.E.  (1959)  Glossary  of  Meteorology.  Amur.  Meteorol.  Soc. , 
Boston.  Mass. 

-It  might  be  argued,  in  view  of  the  power  function  form  of  Kq.  (20),  that  it  would 
bo  preferable  to  assume  a  logarithmic  variation  of  k  with  altitude,  rather 
than  linear.  Practically,  however,  these  different  assumptions  do  not  change 
the  computed  M^  values  significantly  (relative  to  the  inherent  uncertainties  of 
the  overall  analyses). 


by  a  factor  of  two  or  so  within  a  relatively-smal!  descent  distance  of  a  few  thou¬ 
sand  feet.  The  number  density  of  these  new.  larger  particles  also  increased 
downward. 

Such  observations  suggest  that  some  "quantum-type"  of  growth  process  might 
have  been  occurring  in  which  the  largest  snowflakes  were  themselves  colliding  and 
"sticking  together"  to  form  conglomerates  that  were  as  much  as  twice  as  large  as 
those  existent  at  slightly  higher  altitudes.  It  is  also  suspected  that  liquid  water 
droplets  (cloud  droplets)  might  have  been  present  in  the  storm  at  these  altitudes 
which  would  cover  the  falling,  tumbling  snowflakes  with  u  thin  outer  coating  of 
liquid-water,  freezing- water,  or  rime.  Such  coatings  would  greatly  enhance  the 
cohesion  of  these  snowflakes  following  collisions. 

Rased  on  these  observations  and  hypotheses,  large-snow  was  tentatively  de¬ 
fined  to  consist  of  populations  of  snow  particles  having  maximum  particle  sizes 
exceeding  1  mm  equivulent-melted-diameter  (corresponding  to  physical  sizes 
ranging  from  about  2  to  6  mm). 

With  this  defirition,  and  with  that  for  small-snow  previously  stated,  the  tran¬ 
sition  zone  separating  the  storm  regions  of  small-snow  aloft  from  large-snow  be¬ 
low  also  became  defined.  :  Thus,  the  transition  zone  became  the  particular 
growth  layer  in  the  storm  in  which  the  maximum  particle  sizes  of  the  falling  snow 
increased  from  0.5  to  1  mm,  in  equivalent-melted-diameter,  or  from  physical 
sizes  ranging  from  about  1  to  3  mm,  dependent  on  snow  type,  to  sizes  ranging 
from  about  2  to  6  mm. 

The  computations  of  P  and  M  for  the  portions  of  the  missile  trajectories  that 
intersected  the  transition  zone  between  small-snow  and  large-snow  were  per¬ 
formed  in  a  manner  analogous  to  that  described  in  Section  4.5.  3  .  The  precipita¬ 
tion  rate  was  assumed  to  be  linearly  variable  with  altitude  and  to  be  given  by 


'"This  does  not  imply  that  collisions  were  restricted  only  to  the  largest  snowflakes. 
Rather,  it  implies  that  observational  concentration  was  devoted  primarily  to 
these  particular  flakes  within  this  zone. 

:':There  are  pragmatic  reasons,  in  addition  to  any  attempted  physical  justification, 
for  specifying  storm  regions  of  small-snow  and  large-snow,  together  with  a 
transition  zone.  The  size-spectrum  variability  of  snow  in  general  is  simply 
too  large  to  permit  accurate  computations  of  P  and  M  from  empirical  equa¬ 
tions  of  the  Table  2  type  without  some  size  categorization,  such  as  the  two- 
region  categorization  defined  herein.  The  transition  zone  is  necessary  to 
prevent  computational  discontinuities  of  P  and  M  ,  such  as  would  result 
across  the  boundary  between  regions  in  the  absence  of  a  transition  zone. 

What  is  really  needed  to  resolve  certain  of  these;  difficulties  with  snow  is  the 
acquisition  by  snow-type  of  equation  sets  relating  P ,  /, ,  and  M,  which 
properly  reflect  the  functional  dependence  of  these  parameters  on  the  size- 
spectrum  range  of  the  snow  particles.  We  do  not  have  such  equations  at 
present,  nor  the  data  that,  would  permit  Iheir  obtainment. 
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pd  =  PSSB  +  a2d' 


(33) 


where  P.  and  d  are  as  previously  specified,  where  PssD  the  precipitation 
a  k  B 

rate  at  the  base  of  the  small- snow  region  and 


‘  pSSg 


H; 


(34) 


where  P^g^  *s  t^'e  precipitation  rate  at  the  top  of  the  large- snow  region  and  H2 
is  the  total  vertical  depth  of  the  transition  zone. 

The  Mj  values  were  computed  from  the  P^  values  using  the  equation 


(35) 


where  kgg  and  egg  are  the  coefficient  and  exponent  of  Eq.  (20),  that  have  the 
Table  2,  column  2,  ;  alues  pertaining  to  small-snow  where 


^2 


(36) 
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and  where  k^g  and  eLg  are  the  constant  and  exponent  of  Eq.  (20)  having  the 
Table  2,  column  2,  values  corresponding  to  the  type  of  large  snow  observed  in 
the  Wallops  storms. 

4.5.6  LARGE-SNOW  REGION 

The  definition  of  large-snow  was  stated  previously.  The  large-snow  region 
extended  downward  from  the  base  of  the  transition  zone  just  discussed  to  the 
storm  level  at  which  the  falling  snowflakes  first  began  to  melt.  The  base  of  the 
snow  region,  in  other  words,  was  the  altitude  of  the  0°C  isotherm,  or  very  nearly. 

The  P  and  M  values  for  the  missile  trajectory  segments  that  passed  through 
this  region  were  computed  from  the  Table  2  equations  for  the  observed  snow  type. 

4.5.7  MELTING  ZONE 

The  upper  boundary  of  the  melting  zone  was  specified  to  be  the  storm  level  at 
which  the  snowflakes  first  began  to  melt.  The  lower  boundary  was  designated  to 
be  the  level  at  which  all  melting  was  complete,  and  below  which  only  water- 
hydrometeors  (rain)  existed. 


35 


Dependent  on  the  temperature  structure  of  any  given  storm,  it  is  sometimes 
possible  to  have  two  or  more  melting  zones.  Sandwiched  layers  of  above-freezing 
and  below-freezing  temperatures  can  cause  alternate  melting,  refreezing,  and 
remelting  of  the  snow- ice  particles.  Such  multiple  zones  are  not  commonly  ob¬ 
served  but  do  occur  occasionally,  as  in  the  case  of  the  Wallops  storm  of  17  Feb¬ 
ruary  107  2. 

The  melting  zone  is  customarily  detected  by  the  Wallops  Island  radars  (the 
FPS-18,  Spandar,  and  the  X  band).  ''  The  zone,  on  RHI  (range  height  indicator! 
scope  presentations,  appears  as  a  distinctive  "bright  band".  This  bright  band 
occurs  because  the  .‘nowflakes  become  water-coated  during  melting  and  reflect 
microwave  energy  as  if  they  were  large  water  particles  of  snowflake  size.  The 
enhanced  radar  return  is  explained  by  the  fact  that  the  reflectivity  of  liquid  water 
at  microwavelengths  is  substantially  (about  5  times)  larger  than  that  for  ice. 

As  noted  in  R  No.  t,  the  radar  voiume  reflectivity,  rj,  is  measurable  in  the 
melting  zone,  but  the  radar  reflectivity  factor,  Z,  is  indeterminate. 

The  P  and  M  computations  for  the  melting  zone  were  performed  using  equa¬ 
tions  similar  to  those  for  the  other  transition  zones. 

Thus, 


Pd  =  PFSr  +  “id  • 


M, 


rT  " 

H, 


*1 


kr  -  kls 


and 


eR  ’  £LS 


(38) 

(39) 

(40) 

(41) 

(42) 


where  P^  ,  and  d  are  as  previously  defined,  where  RLSg  and  Pr-j,  are  the 
precipitation  rates  at  the  base  of  the  large-snow  region  and  at  the  top  of  the  rain 
region,  respectively,  where  is  the  vertical  depth  of  the  melting  zone,  and 
where  Kj  ^  and  €j  ^  are  the  coefficient  and  exponent  of  Eq.  (20)  that  have  the  par¬ 
ticular  values  listed  in  Table  2,  column  2,  for  the  observed  type  of  large-snow. 

:'The  melting  zone  is  sometimes  not  detected  if  the  storm  has  internal  convective 
structure  consisting  of  cells  of  updraft  and  downdraft  motions. 
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4.  5.  8  MAIN  REGION 


The  rain  region  extended  downward  from  the  base  of  the  melting  zone  to  the 
surface  level.  Rain  was  invariably  present  at  the  surface  level  in  all  Wallops 
storms  to  date  into  which  SAMS  missiles  have  been  fired. 

The  P  and  M  computations  for  the  missile  trajectory  segments  crossing  the 
rain  region  were  handled  differently  than  in  the  other  hydrometeor  regions  dis¬ 
cussed  previously.  The  principal  factor  of  difference  was  that  the  size  distribu¬ 
tion  properties  of  the  rain  at  the  surface  level  were  measured,  for  each  individual 
storm,  by  disdrometer  instruments  located  at  or  near  the  launch  site  of  the  mis¬ 
siles.  Thus,  from  the  measurements  of  P,  M,  and  Z  (described  in  R  No.  1),  the 
particular  P  vs  Z  ,  M  vs  P,  and  M  vs  Z  equations  could  be  established  by  cor¬ 
relation  and  regression  analyses  for  each  storm  for  which  reliable,  quantitative 
disdrometer  data  had  been  acquired. 

The  equations  obtained  for  the  storms  of  the  1971-72  season  of  SAMS  opera¬ 
tions  are  presented  in  Table  2  at  the  bottom.  These  equations  were  used  in  the 

P  and  M  computations  for  the  rain  region  in  the  1971-72  season.  The  equations 
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of  .loss  et  al,  for  widespread  rain  (see  Table  21  were  used  for  the  rain  region 
computations  of  the  1972-73  season. 

4.5.9  SURFACE  LEVEL 

The  precipitation  rate  at  the  surface  level  at  the  launch  time(s>  of  the  mis- 
sile(sl  was  measured  directly  by  means  of  a  tipping-bucket  rain  gauge  located 
near  the  launch  site.  This  was  the  value  of  P  used  at  the  su.face  level. 

The  corresponding  M  value  at  the  surface  level  was  computed  from  Eq.  (201, 
using  the  "disdrometer  values"  of  k  and  e  of  Table  2  for  the  storms  of  the 
197  1-72  seasons,  and  using  the  .loss  values  for  widespread  ”ain  for  the  storms 
of  the  1972-73  season. 

4.6  Trajectory  Computation*  from  the  Surface  Upward 

Because  of  the  need  to  compute  the  integrated  amounts  of  liquid  water  that 
were  intercepted  by  the  missiles  as  a  function  of  their  travel  distance  from  the 
launch  pad  along  the  trajectories,  the  transition  zone  equations  of  Sections  4.5.3, 
4.  5.  5  .  and  4.  5.  7  had  to  be  written  in  modified  form,  such  that  the  P  and  M 
values  were  specified  in  terms  of  distance  measured  upward  from  the  base  of  the 


:  The  equations  obtained  for  the  1970-71  season  are  unknown  to  the  author.  Those 
for  the  1972-73  season  are  available,  but  because  of  the  data  interpretation 
problems  cited  in  R  No.  1,  they  were  obtained  too  late  for  application  to  the 
P  and  \1  computations  for  the  1972-73  season.  These  latter  equations  are 
not  listed  iri  '’’able  2. 
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zones,  rather  than  downward  from  the  top.  The  P  and  M  equations  for  the  first 
transition  zone  above  the  surface,  the  melting  zone,  were  therefore  rewritten  as 


ph  =  pRt  "  a  l  h 


(431 


(44) 


where  h  is  distance  measured  upward  from  the  base  of  the  zone,  where  Oj ,  )3j , 
and  yj  are  as  specified,  respectively,  by  Eqs.  (40),  (41),  and  (42),  and  where  the 
other  parameters  are  as  previously  defined. 

Analogously,  the  equations  for  the  second  and  third  transition  zones  above  the 
surface  were  modified  to 


P 


h 


P 


LS 


a2h 


(45) 


and 


for  *he  small-snow  to  large-snow  transition  zone,  and  to 


Ph  =  PSS  *  a3h 


(46) 


(47) 


and 


Mh  =  (kSS  '  ^3 h)  Ph 


(4SS  "  y3h) 


(48) 


for  the  ice-crystal  to  small-snow  transition  zone,  where  h  ,  as  before,  is  dis¬ 
tance  measured  upward  from  the  base  of  the  zones  and  a2  •  *  ^2  ’  ■  ^3  *  and 

y3  are  specified  by  Eqs. (34),  (36),  (37),  (29),  (31),  and  (32),  respectively.  The 
other  parameters  were  defined  earlier. 


4.7  Integral  of  l,k|uid- Water-Content  Along  the  Missile  Trajectory 

The  equations  described  in  the  preceding  sections  permitted  the  computation 
of  P  and  M  for  the  SAMS  missile  trajectories. 

With  knowledge  of  M  ,  additional  computations  were  performed  to  provide  in¬ 
formation  about  the  ’’trajectory  integral"  of  M  as  a  function  of  the  path  distance, 
R<- ,  from  the  missile  launch  pad.  Thus,  the  values  of 
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(49) 


J  MdRg  = 

0 

were  computed  for  each  radar  (or  interpolated)  data  point  along  the  trajectory  aa 
a  function  of  the  path  distance 


In  these  equations,  ARg^  is  the  incremental  distance  between  the  successive 
data  points,  which  are  numbered  from  i=l  to  i=n  along  the  trajectory,  and  Mj  is 
the  average  liquid-water-content  between  the  successive  points.  The  summation 
is  accomplished  cumulatively,  from  data  point  to  data  point,  up  to  the  storm  top 
level. 


4.8  Profile  Value*  of  the  Hydrometeor  Parameter* 

The  trajectory  values  of  P  and  M  were  tabulated  and  plotted  as  a  function  of 
the  altitude,  z ,  of  the  data  points  above  the  surface  level.  These  tabulations  and 
plots,  with  z  as  the  vertical  coordinate,  will  be  referred  to  herein  as  the  "profile 
values"  or  "profile  plots".  This  method  of  tabulation  and  plotting  facilitates  com¬ 
parisons  of  the  erosion  conditions  encountered  by  the  different  missiles  launched 
into  the  different  storms. 

With  regard  to  the  profile  values  of  the  integral  of  liquid- water- content,  an 
assumption  was  made  regarding  the  relationship  between  the  trajectory  distance 
between  data  points,  ARg  ,  and  the  altitude  distance  (component  distance)  between 
the  same  points,  Az.  These  increments  are  related  as 


ARg  =  esc  6  Az  , 


(51) 


where  4  is  the  elevation  angle  of  the  missile  trajectory  above  the  horizontal  plane 
in  the  vicinity  of  the  data  points. 

This  equation  substituted  in  Eq.  (49)  yields 


1ft.  (esc  ^)j  Az.  , 


(52) 


which  permits  the  integral  on  the  left  to  be  determined  as  a  profile  function  of  the 
altitude 
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z  = 


(53) 


rather  than  as  a  function  of  the  path  distance,  RQ,  as  formerly  with  Eqs.(49)  and 
(50). 

The  trajectories  of  the  Terrier- Recruit  missiles  used  in  the  SAMS  program 
of  tht  1970-73  seasons,  the  portions  between  the  surface  and  storm  top  altitudes, 
can  be  approximated  quite  well  by  a  straight  line  having  an  avers  f,e  elevation 
angle,  4>.  This  permits  Eq.  (52)  to  be  written  as 

Rs 

y  MdRg  ?  esc 
0 


For  the  five  SAMS  missiles  fired  during  the  1971-72  season,  the  average 
elevation  angle  of  the  straight  lines  of  "best-fit  approximation"  to  the  actual  tra¬ 
jectories  was  29.3°.  This  value,  inserted  in  Eq.  (54),  for  yields 


(55) 


MS  l=n 

j*  MdRg  ?  2.04  Az.  . 

°  °  RS 

The  above  equation  was  used  to  compute  the  profile  values  of  \  MdRg  for 


both  the  1971-72  and  1972-73  seasons  of  SAMS  operations,  "'he  values  of  the 
integral  thus  computed  should  not  depart  by  more  than  about  5  percent  at  the 
storm  top  level  of  maximum  value  from  those  computed  from  the  more  rigorous 
Eqs.  (49)  and  (50). 


5.  COMPUTATIONAL  EXAMPLE 

A  typical  example  has  been  selected  that  will  serve  to  summarize  and  illus¬ 
trate  the  computational  procedures  described  in  the  preceding  sections.  The  par¬ 
ticular  example  is  the  missile  flight,  No.Q2-6361,  of  2  February  1973,  launched 
at  1408:30  GMT. 

The  trajectory  values  of  the  radar  integration  signal,  Jj  ,  for  this  flight  are 
shown  in  Figure  1  as  a  function  of  altitude;  see  the  upper  abeissa  scale.  *  The 
solid,  middle  portion  of  the  profile  shows  the  radar-measured  values  of  J_,  for  the 
trajectory  portion  extending  from  just  above  the  ground  clutter  region  near  the 


*The  "bracket  underline"  symbol,  "LJ"  ,  is  used,  as  in  R  No.  1_,  to  specify  the 
decibel-equivalent  value  of  any  given  parameter.  Thus,  ^  =  10  logT. 
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Figure  1.  Profile  of  the  Integration  Signal,  J,.  for 
the  Storm  Trajectory  of  Missile  Q2-6361,  of  2  Feb¬ 
ruary  1973,  The  associated  decibel  values  of  radar 
volume  reflectivity,  j}/  are  indicated  by  the  second 
abscissa  scale.  The  solid,  dashed  and  dotted  por¬ 
tions  of  the  profile  are  explained  in  the  text 


surface  level  to  the  storm  altitude  where  the  value  became  smaller  than  the 
minimum  detectable  value  of  the  radar. 

The  dashed  portion  of  the  profile,  at  the  bottom,  shows  the  ^  values  mea¬ 
sured  along  a  vertical  line  located  about  1.3  km  offshore  of  the  missile  launch  site, 
where  ground  clutter  effects  were  absent  or  minimized,  relative  to  those  at  the 
launch  site.  The  values  were  presumed  to  be  representative  of  the  trajectory 
values. 

The  dotted  portion  of  the  Figure  1  profile,  at  the  top,  shows  the  ^  values  that 
were  assumed  to  apply  to  the  uppermost  part  of  the  missile  trajectory.  The  data 
ooint  at  8.8  km  altitude,  for  which  45.5  db,  was  obtained  by  "special  process¬ 
ing"  of  the  radar  data  recorded  at  the  launch  time  of  the  missile.  The  radar  signal 
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values  for  a  1°  sector  of  the  RHI  scan,  for  an  elevation  angle  of  approximately 
30°,  were  integrated  over  the  sector,  for  each  and  all  pulse-volume  ranges  in  the 
upper  portion  of  the  storm,  to  establish  the  minimum  detectable  value  of  the  in¬ 
tegral  and  the  associate  1  storm  altitude.  *  This  technique,  in  essence,  employs 
integration  to  determine  the  values  of  "weak  signals"  that  are  "obscured  by"  re¬ 
ceiver  noise.  T!  •  data  point  cited  above  was  ascertained  by  this  technique. 

The  storm  top  altitude  on  2  February  1973  was  10.1  km,  as  established  by 
observation  from  the  AFCRL  C-130  aircraft.  The  1  value  at  the  storm  top  was 
assumed  to  be  unity  (that  is,  ^  =  0),  and  the  I  values  between  data  points  in  the 
upper  portion  cf  the  storm  were  assumed  to  be  linearly  variable  with  altitude.  ** 

The  ^  value  at  the  s.mface  level,  specifically  emphasized  in  Figure  1  by  the 
large  black  dot,  was  inferred  from  the  precipitation  rate  measured  by  the  tipping- 
bucket  rain  gauge  at  the  launch  site  at  the  firing  time  of  the  missile,  t 

The  lower  abscissa  scale  in  Figure  1  shows  the  decibel  values  of  the  radar 
volume  reflectivity,^,  that  correspond  to  the  Rvalues,  and  are  related,  as 
specified  by  Eq.  (47'  of  R  No.  1. 

Aircraft  measurements  during  the  storm  of  2  February  1973  established  that 
the  base  and  top  of  the  melting  zone  were  located  at  altitudes  of  about  2.7  and 
3.7  km,  respectively.  This  information  permitted  the  values  of  the  reflectivity 
factor  for  water  hydrometeors,  ZW  *  and  those  for  ice  hydrometeors,  Zj ,  to  be 
computed  from  the  ^  values  of  the  Figure  1  profile,  utilizing  Eqs.(72)  and  (73) 
of  R  No.  1  .  The  resultant  profiles  of  Z^-  and  Zj  for  the  missile  trajectory  are 
shown  in  Figure  2.  The  separate  profiles  have  been  extended  across  the  melting 
zone  itself.  The  values  of  the  reflectivity  factor  are  indeterminant  within  the 
zone,  as  mentioned  earlier.  The  values  should  generally  exceed  those  of  the  Zj 
profile,  but  the  amount  of  the  excess  is  impossible  to  predict  without  making 
theoretical  assumptions  about  the  nature  of  the  melting  processes  in  this  particu¬ 
lar  storm. 

The  aircraft  measurements  and  observations  also  provided  information  about 
the  upper  levels  of  the  storm.  The  large-snow  in  the  region  immediately  above 


•'The  maximum  elevation  angle  of  the  RHI  scans  during  the  stor  i  of  2  February 
1973  was  about  30°. 

**This  assumption  was  used  for  the  1972-73  SAMS  data.  However,  the  assump¬ 
tion  that  J  was  linearly  variable  was  used  for  the  1971-72  data.  It  should  be 
noted  that  if  I  varies  linearly  with  altitude,  does  not,  and  vice  versa.  This 
explains  the  curvature  of  the  Figure  1  profile  between  the  circled  data  points 
above  7  km  altitude. 

tin  other  words,  to  complete  the  profile,  the  J,  value  for  the  surface  level  was 
"back  calculated"  from  H ,  first  utilizing  the  Z  vs  I‘  relationship  of  .loss  et 
al,29  f(>_r  wide«t>'-ead  rain,  as  listed  in  Table  3,  column  1,  and  then  employ¬ 
ing  the  vs  Z  relationship  of  Kq.<72)  of  H  No.  1  . 
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Figure  2.  Profiles  of  the  Radar  Reflectivity 
Factors,  Z\\'  ,  for  Water  Hydrometeors,  and  Zj, 
for  I ce  Hydrometeors,  for  the  Missile  Trajectory 


the  melting  zone  was  composed  primarily  of  "aggregates  of  plates".  The  base  and 
top  of  the  small-snow  to  large-snow  transition  zone  were  located  at  altitudes  of 
about  5.0  and  6.8  km,  respectively.  The  base  and  top  of  the  ice-crystal  to  small- 
snow  transition  zone  were  located  at  about  8.0  and  8.6  km,  respectively.  These 
hydrometeor  regions  and  transition  zones  are  illustrated  in  Figure  3. 

Mso  illustrated  in  Figure  3  is  the  profile  of  precipitation  rate  for  the  missile 
trajectory.  The  P  values  for  the  identified  hydrometeor  regions  and  types  were 
computed  from  the  and  /.[  values  of  the  Figure  2  profiles  using  the  pertinent 
equations  of  Table  2,  column  1.  The  P  values  within  the  transition  zones  were 
assumed  to  be  linearly  variable  with  altitude  as  previously  discussed,  and  as  is 
revealed  in  Figure  3  by  the  straight-line  segments  shown  crossing  the  zones. 
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Figure  3,  The  Hydrometeor  Regions  and  Zones 
Within  the  Storm  and  the  Profile  of  Precipitation 
Rate,  P 


The  M  values  for  the  missile  trajectory  are  illustrated  by  the  solid-line  pro¬ 
file  of  Figure  4.  These  M  values  for  the  different  regions  and  zones  were  com¬ 
puted  from  P  using  the  applicable  equations  of  Table  2,  column  2,  within  regions, 
and  the  interpolation  equations  of  Section  4.6  within  zones. 

The  dashed  curve  of  Figure  4  shows  the  profile  of  J  M  dRg .  The  values  of 
this  integral  were  determined  from  the  T»l  values,  using  Eqs.  (53)  and  (55)  of 
Section  4.8.  The  maximum  value  of  the  integral  for  the  entire  storm  passage 

distance  of  the  missile  is  indicated  by  the  drafted  number  at  the  top  of  the  profile. 

_2 

The  units  of  the  integral  are  gm  m  . 
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Figure  4.  Profiles  of  Liquid -Water- Content,  M, 
and  Integral  of  Liquid-Water,  J M  dRs  ,  for  the 
Missile  Trajectory 


6.  UNCERTAINTIES 

The  erosion  pa.ameters  of  prime  importance  to  the  SAMS  program  at  Wallops 
Island  are  the  profile  values  of  M  and  j M  dRs.  It  follows,  then,  that  questions 
concerning  the  uncertainties  and  error  bounds  of  these  two  parameters  are  also 
of  major  interest. 

The  uncertainties  and  indeterminacies  of  the  radar-measured  values  of  Z ^ 
and  Zj  were  discussed  in  R  No.  1  in  Appendix  B.  These  uncertainties,  for  the 
SAMS  conditions  of  radar  operation  (with  the  FPS-18  radar)  during  the  1971-72 
and  1972-73  seasons,  were  suggested  to  be  of  the  order  of  ±50  percent  to  ±65 
percent,  for  8  range-cell  integration,  dependent  on  the  estimates  made  regarding 
the  uncertainties  of  the  individual  contributing  factors. 
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It  has  been  shown  herein  that  the  M  values  for  the  missile  trajectories  wore 
derived  from  thp  radar-measured  values  of  Z,  either  Zw  or  Zj ,  using  empirical 
equations  of  the  types  listed  in  Table  2.  The  equation  that  specifies  the  direct  de¬ 
pendence  of  M  on  Z  is  Eq.(25),  which,  rewritten,  is 


M 


K 


M 


(25) 


This  equation  shows  that  the  uncertainties  of  M  will  depend  on  the  component 
uncertainties  of  each  of  the  right-hand  factors,  Z,  K^, ,  and  Ejy.*v 

The  total  variation  of  M  resulting  from  component  variation  of  Z  ,  K^,  ,  and 
Em  is  given  by 
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which,  when  the  partial  derivatives  are  evaluated  from  Eq.  (25),  becomes 
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In  ZdE 


(57) 


The  derivatives  dz  .  dK^,  and  dEyj  of  this  equation  may  be  identified  and 
associated  with  the  uncertainty  variations  of  the  parameters  Z,  ,  and  EM. 

The  possible  magnitude  of  these  uncertainties  will  be  discussed  in  the  follow¬ 
ing  sections,  and  estimates  will  be  made  of  the  associated  uncertainties  of  liquid- 
water- content. 

Two  kinds  of  situations  are  postulated  for  uncertainty  estimation.  In  the  first, 
which  will  subsequently  be  referred  to  as  "Situation  l",  it  is  assumed  that  only 
radar-measurement  data  are  available  and  that  auxiliary  surface  and  aircraft  data 
are  not  (except  for  knowledge  of  the  general  category  of  the  hvdrometeors;  that  is, 
whether  they  are  rain,  snow,  or  ice  crystals'.  In  "Situation  2",  in  contrast,  both 
types  ot  data  are  presumed  to  have  been  acquired. 


'  The  fact  that  M  ,  in  the  profile  computations  just  discussed,  was  computed  from 
Eqs. ; 20)  and  (24),  rather  than  from  Eq.  (25),  does  not  negate  this  statement. 
The  equation  sets  of  Table  2  are  self-consistent,  and  the  M  values  computed 
by  either  procedure  will  be  identical  (within  two-place  accuracy). 

**Since  Eq.  (25)  is  an  empirical  equation,  questions  can  be  raised  whether  some 
different  form  of  equation,  other  than  a  power  function,  might  be  more  des¬ 
criptive  of  the  data,  hence  more  accurate.  Such  questions  are  valid,  but 
they  are  beyond  the  scope  of  the  present  report.  The  uncertainty  comments 
of  this  section  are  based  on  the  assumption  that  Eq.  (25)  is  ti  e  best  descrip¬ 
tive  equation  presently  available. 
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6.1  btimated  Uncertainty  of  Liquid-Water-Content 
for  Rain  and  Large-Snow  Under  Situation  I 

In  Situation  1.  with  reference  to  Eq.  (57',  the  dM  uncertainties  of  liquid- 
water-content  will  depend  on  four  component  uncertainties:  (1)  on  the  dKjyj  un¬ 
certainties  of  our  knowledge  of  the  proper  constant  of  Eq.  (25)  that  should  be 
applied  to  the  given  hydrometeor  category;  (2)  on  the  dE^  uncertainties  of  our 
knowledge  of  the  proper  exponent;  (3)  on  the  dZ  uncertainties  of  radar  measure¬ 
ment  that  were  discussed  in  R  No.  1  ;  and  (4)  on  the  dZ  uncertainties  that  reflect 
the  scatter  of  the  Z  values,  used  in  the  regression  analyses,  about  the  regression 
line  obtained  from  the  analysis.  Thus,  with  regard  to  the  latt  'r  two  components, 
the  dZ  parameter  of  Eq.  (57)  is  actually  composed  of  two  parts, 

dZ  =  dZ  +  dZc  .  (58) 

rm  a 

where  dZrm  is  the  uncertainty  of  the  radar- measured  values  of  Z  and  dZg  is  the 
uncertainty,  or  "scatter",  of  the  Z  values  used  in  the  regression  analyses  about 
the  regression  line. 

When  Eq.  (58)  is  substituted  in  Eq.  (57),  we  obtain  the  general  form  of  the  un¬ 
certainty  equation  for  M  ,  that  is, 

dM  =  Km  ZEm  (dZrm  *  dZs)  .  z  JeJ  .  <59, 

In  evaluating  this  equation  for  the  Situation  1  presumption,  it  will  be  advan¬ 
tageous  to  discuss  the  dKM  and  dE^j  component  uncertainties  first.  We  will  then 
review  and  retabulate  certain  of  the  dZrm  information  from  R  No.  t.  Finally,  we 
will  consider  the  dZs  uncertainties  and  present  the  equation-evaluation  results. 

The  ranges  of  the  KM  and  E^j  values  for  rain  and  large-snow  are  shown 
from  the  equations  of  Table  2,  column  3.  In  lieu  of  knowledge  of  the  specific  rain 
or  large-snow  type,  which  is  the  Situation  1  presumption,  it  would  seem  that  we 
would  be  subject  to  KM  and  E variability  at  least  as  large  as  the  ranges  shown 
in  this  table. 

The  Table  2  equations  for  rain  (the  disdrometer  and  Joss  equations)  reveal 
that  Km  varies  from  0.00134  to  0.00624  ,  with  an  arithmetic  mean  value  of 
0.00311  and  a  standard  deviation  of  ±0.  00167  .  The  E^  values  range  from  0.490 
to  0.667;  the  arithmetic  mean  is  0.585;  the  standard  deviation  is  ±0.0516. 

For  large-snow,  the  Table  2  equations  show  that  varies  from  0.004  20 
to  0.00807  ,  with  a  mean  of  0.00616  and  standard  deviation  of  ±0.00140.  The  E]yj 
values  range  from  0.394  to  0.615;  the  mean  is  0.522;  the  standard  deviation  is 
±0.  0807  . 
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It  was  assumed,  for  the  purposes  of  uncertainty  estimation,  that  the  K^j  and 
E|y|  values  of  Eq.  (59),  for  rain  and  large-snow,  would  be  given  by  the  mean  values 
cited  above  and  that  the  dK^  and  dE^  uncertainty  components  of  the  equation 
would  be  specified  by  the  cited  standard  deviation  values. 

Since  only  single  M  vs  Z  equations  are  presently  available  for  the  small-snow 
and  ice-crystal  categories  (see  Table  2),  it  is  impossible  to  judge  the  dK^  and 
dEM  uncertainties  that  might  exist  between  regression  equations  for  the  different 
types  of  small-snow,  or  the  different  types  of  Ice  crystals.  About  all  that  can  be 
said,  at  the  moment,  is  that  the  uncertainties  will  probably  be  as  large,  and  per¬ 
haps  larger,  than  those  for  large-snow. 

We  will  now  consider  the  component  uncertainties,  dZ _ and  dZ„  ,  of 

r  rm  s 

Eq.  (59),  which  involve  the  reflectivity  factor. 

Reference  is  made  to  Table  4,  where  the  Jj,  n  ,  Up  ,  and  Z  values  have 

rl 

been  tabulated  for  the  hydrometeor  categories  and  radar  conditions  described  in 
R  No.  1  (in  Tables  B1  through  B4>.  These  same  radar  conditions  are  assumed  to 
apply  herein,  as  a  starting  point  for  estimating  the  dM  uncertainties.  Thus,  ex¬ 
cept  as  noted  in  the  footnote  below,  the  uncertainty  discussion  of  the  present  re¬ 
port  is  an  extension  of  that  of  R  No.  1.  Likewise,  '  le  uncertainty  tables  (Tables 
4  and  6  herein'  are  extensions  of  Tables  B1  and  B2  of  R  No.  1. 

The  values  of  M  (approximate  values)  listed  in  Table  4  were  obtained  from 
the  tabulated  Z  vrlues  utilizing  the  relations 


M  =  0.0311  Z 


.585 


(60) 


for  rain,  and 

M  =  0.00616  Z-522  ,  (61) 

for  large-snow.  The  constant,  K^j,  and  exponent,  E^j,  values  for  these  equa¬ 
tions  above  are  the  mean  values  cited  previously  in  this  section. 

With  regard  to  the  dZg  uncertainties  of  the  regression  analyses,  concerning 
the  "scatter  of  the  data  values  of  Z  about  the  regression  line”,  reference  is  made 
to  Figures  5,  6,  and  7. 

The  first  two  figures  show  the  M  vs  Z  points  and  regression  lines  for  two 
sets  of  disdrometer  data  (for  rain)  which  were  acquired  during  the  Wallops  storm 
of  22  March  1972.  The  standard  errors  of  estimate  (S.  E.  )  for  the  two  particular 


^Uncertainty  information  about  dZrm  was  presented  in  R  No.  1  for  two  values  of 
range-cell  integration,  nc  =  1  and  nc  =  8,  Just  the  nc  =  8  information  is  con¬ 
sidered  herein,  jince  this  was  the  number  of  range  cells  commonly  included 
in  the  SAMS  trajectory  data  acquired  during  the  1971-72  and  1972-73  seasons. 
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Table  4.  Estimated  Uncertainties  of  Liquid-Water-Content  for  Rain  and  Large-Snow  Under  Situation 
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Figure  5.  M  vs  Z  Data  for  Rain,  Obtained  From  Disdrometer  A  for  the  Wallops 
Storm  of  22  March  1972,  The  regression  line  for  the  data  is  indicated  and  the 
thin,  solid  lines  on  either  side  show  the  bounds  of  the  standard  error  of  estimate. 
The  arrows  show  the  variation  of  Z  about  the  two,  particular  Z  values,  30.92 
and  3092,  which  are  listed  and  used  in  the  uncertainty  computations  of  Tables  4 
and  5 

Z  values  listed  in  Table  4  are  illustrated  in  the  figures  by  the  horizontal  arrows. 
The  presented  data  are  typical  of  disdrometer  data  obtained  for  the  other  storms 
of  the  1971-72  season. 

Figure  7  show's  similar  data  for  large-snow  of  dendritic-stellar  type.  The 

17 

data  are  those  of  Ohtake  and  Henmi,  while  additional  computations  and  the  re¬ 
gression  analysis  were  performed  by  Vardiman.  * 

For  the  uncertainty  estimations  herein,  it  was  assumed  that  the  dZg  values  of 
Kq.  (59),  for  the  case  of  rain,  would  be  given  approximately  by  the  "half  differ¬ 
ence"  of  the  Z  values  indicated  by  the  S.  E.  arrows  of  Figures  5  and  6.  These 
dZs  values  are  listed  in  Tables  4  and  5,  under  the  section  labeled  "second  term 
contributions". 

The  dZs  values  for  large-snow  were  assumed  to  be  given  by  the  "half  differ¬ 
ence"  of  the  Z  values  indicated  by  the  S,  E.  arrows  of  Figure  7.  These  values 
are  also  listed  in  Tables  4  and  5. 


’'Unpublished  work  performed  for  AFCRI.  by  I..  Vardiman  of  Colorado  State  Univer¬ 
sity  (on  reserve  status  from  the  Air  Weather  Service). 
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Figure  6.  M  vs  Z  Data  for  »lain,  Obtained  From  Disdrometer  B  for  the  Wallops 
Storm  of  22  March  1972.  The  regression  line  for  the  data  is  indicated  and  the  thin, 
solid  lines  on  either  side  show  the  bounds  of  the  standard  error  of  estimate.  The 
arrows  show  the  variation  of  Z  about  the  Two,  Particular  Z  Values,  30.92  and 
3092,  which  are  listed  and  used  in  the  uncertainty  computations  of  Tables  4  and  5. 


The  term  contributions  and  total  uncertainties  of  liquid-water-content  are 
shown  in  Table  4  that  were  computed  from  Eq.  (59)  using  the  dKM  ,  dEM  ,  dZrm , 
and  dZg  estimates  stated  above.  The  dM  (total)  uncertainties  listed  in  the  next 
to  the  last  column  were  obtained  by  probable  error  (P.  E.)  summation,  that  is. 


dM  =  yj dM2  +  dM22  +  dM2  +  dM42  .  (62) 

of  the  separate  term  contributions,  dM^,  dMg,  dMg,  and  dM^  ,  of  Eq.  (59);  see 
Eshbach,'*’9  for  example,  for  a  description  of  this  type  of  summation. 

The  Table  4  results  for  the  Situation  i  assumption  reveal  that  the  uncertain¬ 
ties  of  the  liquid-water-content  values  along  the  SAMS  missile  trajectories  range, 
in  percentage  terms,  from  about  *44  percent  to  *111  percent.*  The  greatest 

59.  Eshbach,  O.  W.  (1957)  Handbook  of  Engineering  Fundamentals.  John  Wiley 
and  Sons,  Inc. ,  New  York,  N.  V.  ,  Sixth  Printing 
♦Although  it  is  convenient  to  cite  percentage  uncertainties  for  qualitative  compari¬ 
son  purposes,  it  should  be  noted  that  large  uncertainties,  such  as  the  above 
which  are  expressed  in  percent,  are  highly  non-symmetric  in  their  positive 
vs  their  negative  values.  For  example,  the  positive  value  can  exceed  100 
whereas  the  negative  value  cannot.  The  meaningful  uncertainty  values  of 
Table  4,  also  of  Table  5,  are  the  dM  values  of  the  next  to  the  last  columns. 
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Figure  7.  M  Vs  Z  Data  of  Ohtake*  for  Large-Snow  of  Dendritic-Stellar  Type. 
The  regression  line  for  the  data  is  indicated  and  the  thin,  solid  lines  on  either 
side  show  the  bounds  of  the  standard  error  of  estimate.  The  arrows  show  the 
variation  of  Z  about  the  two,  particular  Z  values,  137.6  and  13,  760,  which  are 
listed  and  used  in  the  uncertainty  computations  of  Tables  4  and  5 


uncertainties  occur  for  large-snow  having  large  liquid-water-content;  the  smallest 
occur  with  rain  having  small  liquid- water-content. 


6.2  Estimated  Uncertainties  of  Liquid-Wster-Contrnt 
for  Rain  and  Large-Snow  Under  Situation  2 


If  quantitative  surface  and  aircraft  measurements  exist  specifying  the  size- 
distribution  properties  of  the  storm  hydrometeors  (and  supplementing  the  radar 
measurements,  which  is  the  Situation  2  presumption),  these  data,  appropriately 
analyzed,  provide  information  about  the  specific  regression  equations  of  \i  vs  Z 
that  pertain  to  each  hydrometeor  type  observed  in  the  storm. 

In  this  situation,  since  there  is  no  doubt  as  to  which  regression  ecuation 
applies,  the  dK^  and  dE^  uncertainty  components  of  Eq.(59)  have  zero  values 
and  the  equation  reduces  to 


dM 


(63) 


♦Unpublished  data,  «ee  footnote  on  p.  50. 
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For  the  same  assumptions  about  dZrm  and  dZg  that  were  stated  previously, 
the  term  contribution  and  total  uncertainties  of  liquid-water-content  are  shown  in 
Table  5  which  were  computed  from  Eq.  (60).  Probable-error  summation  of  the 
term  contributions  was  used  to  obtain  the  dM  values,  as  in  the  case  of  the  Table 4 
computations. 

A  comparison  of  the  dM  values  of  Tables  4  and  5  reveals  that,  when  dis- 
drometer  data  are  available  to  supplement  the  SAMS  radar  measurements,  the 
uncertainties  of  the  trajectory  values  of  liquid-water- content  in  the  rain  region  of 
the  storm  are  reduced  by  about  ±4  percent  to  ±16  percent,  relative  to  the  uncer¬ 
tainties  of  the  radar  measurements  alone.  The  accuracy  of  measurement,  in 
other  words,  is  improved  by  these  amounts.  Similarly,  when  aircraft  data  are 
acquired  to  supplement  the  radar  measurements,  the  accuracy  of  the  liquid-water- 
content  values  in  the  large- snow  region  of  the  storms  is  improved  by  some 
±  13  percent  to  ±34  percent,  relative  to  the  uncertainties  of  the  radar  measure¬ 
ments  alone.  The  accuracy  enhancement  to  be  anticipated  from  the  acquisition  of 
aircraft  data  in  the  small-snow  and  ice-crystal  regions  of  the  storm  is  probably 
as  large,  or  larger,  than  these  last  cited  values. 

6.3  Comments 

Two  comments  are  pertinent  concerning  the  Table  4  and  Table  5  results. 

It  must  be  emphasized,  first,  as  was  also  emphasized  in  R  No.  1,  that  the 
table  values  are  merely,  first,  gross  estimates  of  the  possible  uncertainty  values. 
Additional,  better  data  must  be  acquired  and  more-careful  analyses  performed  to 
establish  the  definitive  values  of  the  uncertainties. 

Second,  it  should  be  noted  that  an  a-priori  assumption  was  made,  under  the 
Situation  1  hypothesis  (of  radar  measurements  only),  that  knowledge  existed  con¬ 
cerning  the  general  category  of  hydrometeors  that  were  producing  radar  echo, 
whether  they  were  rain,  large-snow,  small-snow,  or  ice  crystals.  In  actual  fact, 
it  is  impossible,  by  radar  means  alone,  to  differentiate  the  latter  three  categories 
of  ice  hydrometeors  which  may  be  present  at  storm  altitudes  above  the  melting 
level.  Aircraft  measurements  are  required  to  do  this.  The  point  here,  which  is 
an  important  point  that  could  not  be  verbalized  previously,  is  that  the  true  uncer¬ 
tainties  of  the  liquid-water-content  values,  as  assessed  from  radar  measure¬ 
ments  only,  for  storm  altitudes  above  the  melting  level,  may  be  substantially 
larger  than  the  single  category,  large-snow  uncertainties  shown  in  Table  4.  The 
dKjyj  and  dE^  component-uncertainties  that  enter  Eq.  (59)  will  actually  be  larger 
than  the  Table  4  estimates  when  aircraft  information  is  lacking,  because  even  the 
hydrometeor  categories  themselves  cannot  then  be  identified  (the  dK^  and  dE^ 
uncertainties  would  reflect  the  larger  "category  differences"  of  Table  2,  rather 
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than  the  "type  differences"  within  the  category).  The  uncertainties  of  liquid-wator- 
content  would  be  correspondingly  increased,  in  such  case,  conceivably  by  as  much 
as  a  factor  of  two,  dependent  upon  the  particular  storm  conditions. 

These  larger,  possible  uncertainties  are  mentioned  because  some  of  the 
SAMS  storm  data  obtained  during  the  1970-73  seasons  consisted  essentially  of 
"radar  measurements  only".  This  would  include  situations  in  which  aircraft  data 
could  not  be  acquired  that  were  representative  of  the  hydrometeor  conditions 
along  the  missile  trajectory  at  firing  time  as,  for  example,  in  the  case  of  missiles 
launched  into  convective-type  storms  having  appreciable  spatial-temporal  vari¬ 
ability. 

6.4  CnrerUintiea  of  the  Integral  of  Liquid-Water-Content 

The  dM  uncertainties  discussed  above  are  the  uncertainties  that  pertain  to 
single-point  measurements  along  the  missile  trajectory  which  have  been  spatially 
averaged  over  the  volumetric  extent  of  the  radar-integration-volume  and  time 
averaged  throughout  the  duration  of  the  integration  period  [see  Eqs.(A3)  and  (A4) 
of  R  No.  lj.  Additional  spatial  averaging  is  performed  when  the  radar  data  are 
processed  by  computer,  as  was  the  case  with  the  SAMS  data  of  the  197  2-73  season. 

The  cloud-physics  parameter  of  primary  concern  to  the  SAMS  ABRES  pro¬ 
gram  is  the  integral  of  liquid-water- content  along  segment  portions  of  the  mis¬ 
sile  trajectory  which  correspond,  in  length  extent,  to  the  resolution  distance  of 
the  erosion  information  that  is  telemetered  to  the  ground  from  the  sensors  aboard 
the  missile.  Detailed  correlations  of  erosion  parameters  vs  cloud  physics  param¬ 
eters  must  necessarily  be  concerned  with  a  "matching  corresnondence"  of  the 
"length  resolution"  of  the  separate  parameters. 

Thus,  it  is  the  path  integral  of  liquid-water-content  at  particular  resolution 
distances  which  is  the  parameter  of  final  cloud-physics  interest  to  SAMS;  and 
the  "accuracy  level"  of  the  AFCRL  measurements  performed  for  SAMS  is  directly 
related  to  the  uncertainties  involving  this  integral.  These  uncertainties  should 
never  exceed  the  value-estimates  discussed  in  the  preceding  sections,  which  per¬ 
tain  to  single-point  measurements  of  M.  Rather,  the  uncertainties  will  be  smaller, 
in  degree  dependent  on  the  nature  and  details  of  such  things  as  the  radar  integra¬ 
tion,  the  computer  averaging,  the  erosion  instrumentation  of  the  missile,  and  the 
analysis  methods  that  were  employed  in  the  particular  Wallops  storms  of  prior 
SAMS  investigation.  The  geometry  of  the  radar  scan  pattern  and  sampling  volume 
relative  to  that  of  the  missile  trajectory  also  enters  the  problem. 

These  matters  should  be  investigated,  if  definitive  information  regarding  un¬ 
certainties  is  desired.  It  is  beyond  the  scope  of  the  present  report,  however,  to 
attempt  such  consideration. 
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7.  SIZE-DISTRIBUTION  INFORMATION  FOR  THE  DIFFERENT  HYDROMETEOR  TYPES 


The  nature  and  amount  of  the  nose-tip  erosion  experienced  by  the  SAMS  mis¬ 
siles  are  influenced  by  the  types  and  size-distribution  properties  of  the  hydro- 
meteors  along  the  missile  trajectories.  This  information,  as  mentioned  previ¬ 
ously,  is  acquired  by  aircraft  measurements. 

Since  the  aircraft  measurements  are  made  subsequent  to  the  launch  times  of 
the  missiles  and  in  spatially  different  parts  of  the  storms  than  those  penetrated  by 
the  missiles,  the  size-distribution  information  obtained  is  not  specifically  appli¬ 
cable  to  the  missile  trajectories.  It  is  merely  representative  of  the  general  storm 
conditions  during  the  time  period  following  the  missile  firing. 

The  size-distribution  analyses  that  have  been  performed  on  the  aircraft  data 
for  the  SAMS  storms  to  date  have  been  relatively  few.  There  were  various  in¬ 
strumentation  failures  during  the  aircraft  flights,  with  consequent  loss  of  data. 

This  limited  the  available  samples.  Also,  the  analyses  are  painstaking  and  time- 
consuming;  hence  the  samples  must  be  selected  with  care  and  be  limited  in  num¬ 
ber.  (These  results  will  be  presented  in  the  forthcoming  data  reports  of  the  series.) 

In  view  of  these  problems  of  the  non-specificity  and  limited  number  of  the  air¬ 
craft  samples  for  the  past  storms  and  for  purposes  of  erosion  assessment,  there 
is  a  need  for  a  method  of  estimating  the  approximate  size  distribution  properties 
of  the  hydrometeors  along  the  missile  trajectories.  One  such  method  will  be  des¬ 
cribed  in  this  section  which  concerns  rain  drops,  snow  aggregates,  and  ice  crys¬ 
tals  of  precipitation  size  (> 80  microns,  approximately).  Methods  for  estimating 
the  size  spectra  of  water  droplets  of  cloud  size  F'80  microns)  will  be  discussed 
later  in  future  reports. 

5  2  4 

Laws  and  Parsons,'  Wexler,  and  Marshall  and  Palmer,  as  mentioned  in 

Section  3.1  ,  demonstrated  that  the  variation  of  the  number  concentration  of  rain¬ 
drops  with  diameter,  D,  could  be  specified  approximately  by  a  distribution  func¬ 
tion  of  exponential  type, 

N  =  N  e‘AU  .  (64) 

o 

3  4 

From  the  data  of  Marshall  et  ai,  Marshall  and  Palmer  demonstrated  that 
-3  -1 

N  had  the  value  8000  (m  mm  )  and  that  A  was  functionally  dependent  on  the 

0  - 1 
precipitation  rate,  P  (in  mm  hr  ),  as 

*  A  1  1,-0. 21  -1 

A  =  4 . 1  P  mm.  (65) 

14 

Gunn  and  Mr"-shall  used  this  same  distribution  function  to  describe  the 
size-distribution  properties  of  aggregate  snow  particles.  The  diameter  D  was 
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replaced  by  the  equivalent-melted-diameter  and  they  ascertained  from  their  data 
analyses  that 

Nq  =  3.8  X  103  P"-87  m-3  mm'1  .  (66) 

and 

A  =  2.55  P"*48  mm"1  .  (67) 

Other  investigators  have  used  the  distribution  function  of  Eq.  (64)  in  conjunc¬ 
tion  with  various  data-dependent  assumptions  about  the  equations  for  N0  and  A  . 
Certain  of  these  investigators  are  identified  in  Table  6  and  their  particular  equa¬ 
tions  are  listed. 


Table  6.  The  Equations  (Data  Supported)  of  Various  Investigators  Regarding 
the  N0  and  A  Terms  of  the  Exponential  Distribution  Function,  Eq.  (64)  Herein 


Equations  for 

Investigators 

Precipitation 

Category 

No 

-3  -1 

m  mm 

A 

mm’1 

4 

Marshall  and  Palmer 

Continuous  Rain 

Nn  =  8,000 

A  —  4  i  p~ 

Joss  et  al40 

Drizzle 

N0  =  30,000 

A  =  5.7  P~’21 

Widespread,  Con¬ 
tinuous  Rain 

N0  =  7,000 

A  =  4.1  P-’21 

Thunderstorm, 

Rain 

Nq=  1,400 

A  =  3.0  P"*21 

30 

Sekhon  and  Srivastava 

Thunderstorm, 

Rain 

N0  =  7,000  P'37 

A  =  3.8  P"*14 

Ramana  Murty  and 

Gupta  (1959) 

Monsoon  Rain, 
Orographic 

N  =  7,500  P"'25 

0 

A  =  0.  81  p-27 

Monsoon  Rain, 
Non-Orographic 

21 

NQ=  1.500  P#tfJ 

17 

A  =  1.14  P#lf 

14 

Gunn  and  Marshall 

Aggregate  Snow 

N  =  3,800  P",8? 

0 

A  =  2.55  P"-48 

23 

Sekhon  and  Srivastava 

Aggregate  Snow 

Nq  =  2,500  P"-94 

A  =  2,29  P"‘45 

In  the  development  herein,  we  will  depart  somewhat  from  the  work  of  these 
previous  investigators.  We  will  use  the  distribution  function  of  Eq.  (64),  but  will 
obtain  expressions  for  N0  and  A  that  are  written  in  terms  of  the  liquid-water- 
content,  rather  than  P.  In  essence,  we  will  follow  the  development  of  Sekhon 
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30 

and  Srivastava.  Additionally,  however,  we  will  develop  equations  which  are 
consistent  with  the  M  vs  Z  relationships  of  Table  2. 

The  liquid-water-content  of  the  drops  or  particles  specified  by  Eq.  1641  for 
all  diameters  (or  equivalent  melted  diameters!  of  the  drops  or  particles,  from 
D  =  0  to  D  =  «> ,  is  given  by 

■y.  00 

M  =  |  X  10‘3  pw  J  N  D3  dD  =  |  X  10‘3  pw  No  J  e'AD  D3  dD  ,  (68) 


-3  -3 

where  pw  is  the  density  of  water,  in  gm  cm  ,  and  M  is  given  in  gm  m  .  On 
integration,  this  equation  becomes 


M  =  |  X  10 


3  <VNorMI 


169) 


where  I"(4)  is  the  gamma  function  of  4. 

The  radar  reflectivity-factor  for  the  distribution  of  drops  or  particles  des¬ 
cribed  by  Eq.  (64)  is,  for  all  diameters  from  D=  0  to  D  =  »  , 


OO  00 

Z  =  f  ND6dD  =  N  f  e'AD  D6  dD  ,  (70) 

Jo  °  Jo 


which,  integrated,  yields 


Z 


No  r(7)  6  -3 

- - —  mm  m 

A? 


(71) 


where  T(7)  is  the  gamma  function  of  7  . 

The  parameter  A  can  be  eliminated  between  Eqs.(69)  and  (71)  to  obtain 


N 


o 


6  X  103  M 
*r<4)  pw 


6  X  103  r(7)  M 

ff  pw  n 4)  z 


(72) 


which,  from  knowledge  that  pw  =  1.0,  T(4)  =  6,  and  T(7)  =  720,  reduces  to 


N 


o 


4.46  X  109 


-3 


m 


-1 

mm 


(73) 


Alternately,  Ny  can  be  eliminated  between  Eqs.  (69)  and  (71)  with  the  result, 
after  evaluation  of  pw,  T(4)  and  T(7), 

*In  practice,  for  practical  equation  accuracy,  integration  to  D  =  «  implies  inte¬ 
gration  to  some  maximum  diameter  approximately  ^  15/A. 
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A 


(74) 


=  61.2 


m 


>1/3 


mm 


-1 


The  empirical  equations  relating  M  and  Z  were  discussed  in  Section  4.2 . 
It  was  pointed  out  that  the  equations  are  of  the  form 


M 


(75) 


and  that  the  particular  values  of  KM  and  E^  for  each  hydrometeor  category  (rain, 
large-snow,  small-snow,  and  ice  crystals)  and  type  were  obtained  from  regres¬ 
sion  analyses  performed  on  data  pertaining  to  the  individual  types.  The  values 
used  in  the  SAMS  program  were  listed  in  Table  2. 

When  Eq.  (75)  is  substituted  into  Eq.  (73)  first,  and  into  Eq.  (74)  second,  wc 
obtain 


4 

9  3  E.. 

N  =  4.46  X  10a  K„  M 
o  M 


7EM-4 


Te 


M 


M 


-3  m  -1 
m  mm 


and 


Em'1 


3  E 


A  =  61.2  K 


M 


Te 


M 


M 


M 


mm 


(76) 


(77) 


which  permit  the  specification  of  the  NQ  and  A  values  for  each  hydrometeor  type 
(from  knowledge  of  KM  and  E^j)  as  functions  of  the  liquid-water-content.  These 
equations,  for  the  hydrometeor  types  defined  for  the  SAMS  program,  are  listed  in 
Table  7.  * 

The  Nq  and  A  equations  of  this  table,  when  substituted  into  the  distribution 
function  of  Eq.  (64),  provide  approximate  information  about  how  the  number  con¬ 
centration  of  the  hydrometeors  varies  with  the  drop  or  particle  diameter  (equiva¬ 
lent- melted-diameter)  for  any  given  value  of  liquid-water-content.  Information  is 
also  provided  about  the  distribution  of  liquid-water-content  within  the  population, 
by  diameter  size.  This  distribution  function,  which  is  the  third  moment  of  the 
function  for  number  concentration,  is 

Md  =  f  X  10*3  pw  Nq  e'AD  D3  gm  m'3  mm"1  .  (78) 

♦The  equations  are  also  shown  in  Table  8,  written  in  terms  of  the  precipitation 
rate,  P  (accomplished  by  use  of  the  P  vs  M  relationships  of  Table  3).  These 
transformed  equations  permit  direct  comparisons  with  the  equations  of  the 
other  investigators  shown  in  Table  6. 
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Table  7.  The  Equations  for  N0 ,  A,  and  Nj  for  the  AFCRL  Hydrometeor  Types  Defined  for  SAMS, 
Written  in  Terms  of  the  Liquid-Water-Content  (See  Table  2  and  Eqs.  (761.  (77),  and  (80)  of  the  text 
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♦This  type  of  snow  was  defined  and  used  only  during  the  1971-72  SAMS  season. 


Table  8.  The  Equations  for  N0  ,  A,  and  N-j-  for  the  AI'CRL  Hydrometeor  Types  Defined  for  SAMS,  Written 
in  Terms  of  the  Precipitation  Rate,  Which  Permits  Comparisons  With  the  Equations  of  Table  6.  The  equa¬ 
tions  listed  were  derived  from  the  Table  7  equations  using  the  M  vs  P  relations  listed  in  Table  2 


6 


*This  type  of  snow  was  defined  and  used  only  during  the  1971-72  season. 


The  N.  and  A  terms  of  this  equation  are  likewise  specified,  for  each  hydrometeor 
type,  by  reference  to  Table  7.  The  independent  variable  of  the  above  equation  is 
D;  the  other  parameters  on  the  right  are  constants  for  any  given  value  of  total 
liquid- water- content. 

1'he  values  of  the  radar  reflectivity  factor  for  the  population  (the  values  per 
diameter  bandwidth'  are  distributed  according  to  the  sixth  moment  of  the  function 
for  number  concentration,  or  as 

Zn  =  N  e‘AD  D6  mm6  m'3  mm-1  .  (79) 

1j  o 


Other  characteristic  or  useful  parameters  of  the  distributions  can  be  specified 
analytically.  For  example,  the  total  number  concentration  of  the  drops  or  particles, 
of  all  sizes  within  the  populations,  is  given  by 


2  E 


M 


Nt  =  7.29  X  10 


K 


\\1 


M 


'M 


m 


-3 


(80) 


(from  integration  of  Fq.  (64)  and  substitutions  from  Eqs.(76)  and  (77)].  The 
and  values  for  the  particular  hydrometeor  tvpes  are  listed  in  Table  7  (also 
in  Table  2). 

The  "median  volume  diameter"  of  the  distribution  function  for  liquid-water- 
content  is  a  parameter  of  common,  conventional  reference.  This  diameter,  first 
derived  by  Atlas60  is, 

Dq  =  3.67/A  mm  .  (81) 

The  "modal  diameters"  of  the  Vlj-j  and  distributions  are  also  character¬ 
istic  parameters  of  interest.  These  diameters,  which  specify  the  peak-value 
points  of  the  bandwidth  values  of  liquid-water- content  and  reflectivity  factor  are 
given  respectively  by 

D'  =  3/A  mm  (82) 

and 

D*  =  6/ A  mm  .  (83) 

z 

These  expressions  were  derived  by  differentiating  Fqs.(78)and  (79)  and  de¬ 
termining  the  diameter  values  corresponding  to  maximum  and  Zp. 

60.  Atlas,  I).  (1953)  Optical  extinction  by  rainfall,  .1.  Mcteorol.  10:486-488. 
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Each  of  the  above  Eqs.  (64)  and  (78)  through  (82),  can  be  evaluated  for  any 
given  point  along  the  SAMS  missile  trajectories.  The  hydrometeor  regions  of  the 
storms  and  the  hydrometeor  types  within  regions  are  established  by  aircraft  ob¬ 
servation,  as  discussed  earlier  and  as  is  illustrated  in  Figures  3  and  4.  This 
information,  by  reference  to  Table  7,  permits  the  determination  of  the  NQ  and  A 
equations  for  the  given  hydrometeor  type,  which  in  turn  permits  the  specification 
of  the  distribution  functions  of  Eqs.  (64),  (78),  and  (79).  These  equations,  also 
Eqs.  (80)  through  (83),  can  then  be  evaluated  for  the  particular  liquid-water-content 
value  that  pertains  to  the  selected  trajectory  point. 

A  specific  example  may  be  helpful. 

Consider  the  profile  plot  of  liquid-water-content  of  Figure  4,  for  the  SAMS 

missile  flight  of  2  February  1973.  Consider  a  point  on  this  profile  which  lies 

within  the  "large-snow  region"  of  the  storm  and  is  located  at  a  trajectory  altitude 

*3 

of  3.8  km.  The  liquid-water-content  value  at  this  altitude  is  0.3  gm  m  (by  con¬ 
venient  choice). 

The  large-snow  on  this  day  was  determined  to  be  of  the  type  called  "aggre¬ 
gates  of  plates".  The  NQ  and  A  equations  for  this  hydrometeor  type  are  as 
shown  in  Table  7.  These  equations,  when  substituted  in  Eqs.  (64),  (78),  and  (79) 
provide  the  particular  distribution  functions, 

N  =  27,700  e'4'12  D  m"3  mm-1  ,  (84) 

for  number  concentration, 

Md  =  14.5  e'4,12  D  D3  gm  m*3  mm”*  ,  (85) 

for  liquid-water-content,  and 

ZD  =  27,700  e'4,12D  D6  mm6  m"3  mm'1  ,  (86) 

for  the  radar  reflectivity  factor. 

Plots  of  these  equations  are  presented  in  Figure  8,  to  illustrate  their 
characteristics.  The  upper  plot  shows  the  size  distribution  of  the  number 
concentration  of  the  snow  aggregates.  The  total  number  of  the  aggregates 
within  the  population,  N-j- ,  is  indicated.  The  plot  in  the  middle  shows  the 
distribution  of  the  liquid-water-content;  the  plot  at  the  bottom  shows  the  dis¬ 
tribution  of  the  radar  reflectivity  factor.  The  integrated  total  value  of  the 

6  “3 

radar  reflectivity  factor  for  the  trajectory  point  was  990  mm  m  ,  which  value 
is  consistent  with  the  M  vs  Z  equation  of  Table  2,  for  large-snow,  of  type  LSg . 
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Figure  8.  Plots  IllusLrating  the  Type  of  Size- 
Spectrum  Information  Ttat  Can  Be  Obtained 
From  the  Non-Truncatec  Distribution  Func¬ 
tions  (Discussed  in  the  Text)  When  They  Are 
Applied  to  Any  Given  Altitude  Point  Along  the 
SAMS  Missile  Trajectories,  Excluding  Points 
Located  Within  the  Melting  Zone 


The  median-volume-diameter  of  the  population,  D  [see  Eq.  (81)],  is  indicated  on 
the  middle  plot;  the  modal  diameters,  D'  and  D'  |see  Eqs.(82)  and  (83)],  are 
indicated  on  the  middle  and  bottom  plots.  Two  abscissa  scales  are  shown.  The 
upper  one  indicates  the  equivalent  melted  diameters  of  the  snow  aggregates;  the 
lower  one  indicates  the  approximate  average  physical  dimensions  of  the  particles 
[see  Eq.  (89)  for  the  scale  relationship  assumed]. 

The  total  number  of  snow  particles  of  all  sizes  in  this  population  example  is 
N.j,  =  6720  from  Eq.  (80)  and  the  and  E^  values  of  Table  7.  The  median  and 
modal  diameters  of  toe  Mp  distribution  are  Do  =  0.89  mm  (equivalent  melted 
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diameter'  and  D' =  0.73  mm  from  Eqs,  (81)  and  (82).  The  modal  diameter  of  the 
ZD  distribution  is  Dz'  =  1.46  mm  from  Eq.(83).  These  diameters  are  indicated 
in  the  Figure  8  diagrams. 

This  example  demonstrates  the  type  of  size-distribution  information  that  can 
be  established  for  any  trajectory  point  of  the  Figure  4  profile  that  is  located 
within  one  of  the  four  (possible)  hydrometeor  regions  specified  herein.  Similar 
information  can  also  be  established  for  two  of  the  three  transition  zones,  the  ex¬ 
ception  being  the  melting  zone.  The  and  Ej^  values  within  the  zones  above 
the  melting  level  can  be  interpolated  with  altitude  above  the  zone  bases,  as  dis¬ 
cussed  in  Sections  4.5  and  4.6 .  The  N0  and  A  equations  for  any  given  height 
within  the  zones  can  then  be  determined  by  direct  evaluation  of  Eqs.  (76)  and  (77), 
which  permits  the  specification  of  the  distribution  functions  of  Eqs.  (64),  (78),  and 
(79);  also  of  the  auxiliary  Eqs.  (80)  through  (83). 

Distribution  function  and  totals  information  is  provided  in  Figures  9  through 
17  for  all  of  the  hydrometeor  types  observed  in  the  storms  of  the  1971-73  SAMS 
seasons.  The  upper  diagrams  show  the  distribution  of  the  number  concentration 
of  the  drops  or  particles  for  three  different  values  of  liquid-water-content,  which 
are  indicated.  The  middle  diagrams  show  the  distribution  of  liquid-water-content 
with  the  equivalent-melted-diameter  of  the  drops  or  particles,  and  the  bottom 
diagrams  show  the  distribution  of  the  radar  reflectivity  factor.  The  total  numbers 
of  the  drops  or  particles  in  the  distributions,  NT ,  are  listed,  as  are  the  total 
Z  values.  The  median-volume-diameter  points  of  the  distributions  (the  DQ  points, 
see  Eq.  (81)]  are  shown  by  the  large  black  dots.  Two  abscissa  scales  have  been 
drawn  for  the  figures  pertaining  to  ice  hydrometeors.  The  upper  scale  shows  the 
equivalent  melted  diameters  of  the  particles;  the  lower  scale  indicates  the  approxi¬ 
mate  physical  dimensions  (average  dimensions)  of  the  particles.  The  scale  equa¬ 
tions  relating  the  physical  length,  l,  to  the  equivalent  melted  diameter,  De,  are 
the  following; 

De  =  0.44  f,0'84  .  (87) 

for  ice  crystals,  type  Cj,  (see  Table  2  for  type  identification) 

De  =  0.324  t0-805  ,  (88) 

for  small-snow,  type  SSg ,  and 


D 


e 


0.40 


j  0.875 


(89) 
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Figure  9.  Distribution-Function  Plots  for  Type  Cj 
Ice-Ci  -stals  for  Three  Different  Values  of  Liquid- 
Water-Content.  The  plots  of  the  upper  diagram 
show  the  size  distribution  of  the  number  concen¬ 
tration  of  the  ice-crystals  for  liquid-water-content 
values  of  0.01  ,  0.05  ,  and  0.1  gm  m"3  .  The  num¬ 
ber  total,  Nt  ,  of  the  ice  crystals  within  the  popula¬ 
tions  is  also  indicated.  The  plots  of  the  middle 
diagram  show  the  distribution  of  the  liquid- water- 
content  by  crystal  size  for  the  same  values  of  total 
liquid-water-content  as  cited  above.  The  median- 
volume-diameters,  D0 ,  are  indicated  by  the  large 
black  dots.  These  diameter  points  are  also  Lhown 
in  the  other  distributions  as  well.  The  plots  of  the 
bottom  diagram  show  the  size  spectra  of  the  radar 
reflectivity  factor.  The  total  integrated  values  of 
the  radar  reflectivity  factor  arc  additionally  shown. 
Two  abscissa  scales  have  been  drafted  on  the  fig¬ 
ure.  The  upper  scale  shows  the  equivalent  melted 
diameters  of  the  ice-crystals;  the  lower  scale  pro¬ 
vides  approximate  information  about  the  average 
physical  sizes  of  the  crystals.  The  scale  relation¬ 
ship  is  specified  by  Kq.  (87) 
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for  large-snow,  type  LSg ,  The  latter  equation  pertains  strictly  to  the  large-snow 
type  LSg ,  but  was  used  for  the  type  LS^  of  Figure  11,  and  the  type  LSj, ,  of  Fig¬ 
ure  12,  in  lieu  of  other  knowledge.  The  equations  were  developed  by  Cunningham 
from  literature  information  and  empirical  data. 

The  distribution  plots  of  Figures  9  through  17  will  serve  to  inform  the  reader 
of  the  general  nature  of  the  size  distribution  characteristics  and  differences  that 
exist  among  the  different  hydrometeor  types  that  have  been  discussed  in  ihis  report. 
Several  comments  are  pertinent  before  concluding  this  section, 
it  should  be  noted' that  the  size-distribution  information  provided  by  the  tech¬ 
niques  and  equations  described  herein  is  only  a  first,  gross  approximation. 
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Figure  10.  Distribution-Function  Plots  for 
Small-Snow,  Type  SS<-  ,  for  Three  Different 
Values  of  Liquid- Water-Content.  The  plots  are 
similar  to  those  described  in  the  Figure  9  cap¬ 
tion.  Note  that  the  scale  of  equivalent  meltec! 
diameter  extends  to  2  mm,  which  differs  from 
the  Figure  9  scale 
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There  are  two  major  reasons  for  this.  First,  the  actual  hydrometeors  in  the 
Wallops  storms  are  not  necessarily  size-distributed  in  an  exponential  manner. 
Observed  distributions  reveal  substantial  individual-sample-departures  from 
exponentiality,  although,  on  the  average,  any  large  number  of  samples  will  re¬ 
veal  distinct  tendencies  toward  exponentiality.  Second,  in  the  development  herein, 
which  follows  the  more  or  less  conventional  treatments  of  the  past,  the  distribu¬ 
tion  function  of  Eq.  (64)  is  integrated  between  D  =  0  and  D  =  «.  *  This  is  physi¬ 
cally  unrealistic,  of  course,  since  drops  or  particles  of  zero  size  do  not  exist. 
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Figure  11.  Distribution- Function  Plots  for 
Large-Snow,  Type  LS3  ,  for  Three  Different 
Values  of  Liquid-Water-Content.  The  plots  are 
similar  to  those  described  in  the  Figure  9  cap¬ 
tion.  Note  that  the  liquid-water-content  values 
are  0.1,  0.2,  and  0.3;  also  that  the  scale  of 
equivalent  melted  diameter  extends  to  5  mm 


*It  can  be  shown  that,  for  practical  purposes,  integration  to  D  =  »  is  equivalent 
to  integrating  to  n  maximum  diameter  defined  by  Dm  t  15/A  ,  or  [see  Eq.  (82)] 

Dm  *  5  D'. 
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and  since  there  is  also  an  upper  diameter  limit  of  the  sizes  of  the  largest  drops 
or  particles.  The  distributions  are  obviously  truncated  at  some  minimum  diam¬ 
eter  value,  also  at  some  maximum  diameter  value.  In  addition  to  this  physical 
truncation  of  the  hydrometeor  populations,  truncation  can  also  occur  instrumen- 
tally,  since  any  given  aircraft  instrument  may  be  incapable  of  sensing  drops  or 
particles  smaller  than  a  particular  minimum  size,  or  larger  than  some  particular 
maximum  size.  These  truncation  problems,  which  are  of  vital  concern  to  SAMS, 
are  currently  being  considered  at  AFCRL.  A  double  truncation  model  has  been 

developed  which  extends  and  simplifies  certain  of  the  equations  of  Sekhon  and 
30 

Srivastava  and  particularizes  them  for  SAMS  application.  This  model  will  be 
described  and  illustrated  later  in  the  report  series. 

The  following  statements  may  be  made  about  the  physical  deficiencies  of  the 
non-truncated  equations  presented  in  this  section  and  illustrated  in  Figures  8 
through  17.  For  any  given  hydrometeor  type  and  liquid  water  content,  the 
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Figure  12.  Distribution-Function  Plots  for 
Large-Snow,  Type  LS^ ,  for  Three  Different 
Values  of  Liquid-Water-Content.  The  plot  scal¬ 
ing  and  liquid- water- content  values  plotted  are 
identical  to  Figure  11 
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non-truncated  equations  (compared  *o  the  more  realistic  truncated  equations)  will 
tend  to  (1)  overestimate  the  values  of  A,  (2)  overestimate  the  values  of  NQ, 

(3)  overestimate  the  number  of  the  smallest  drops  or  particles,  (4)  underestimate 
the  number  of  the  largest  drops  or  particles,  and  (5)  overestimate  the  total  num¬ 
ber  of  the  drops  or  particles.  '' 

Size  distribution  information  about  the  hydrometeors  within  the  melting  zone 
is  not  provided  in  this  report  but  will  be  furnished  later.  The  hydrometeors  in 
this  zone  are  composed  of  fully-melted  water  drops  and  water-coated  ice  particles 
in  varying  proportions,  dependent  on  the  temperature  structure  of  the  atmosphere 
and  the  altitude  within  the  zone.  Spectrum  assessments  for  both  types  of  hydro¬ 
meteors  will  be  presented  as  functions  of  altitude,  decreasing  downward  from  the 
upper  boundary  of  the  melting  zone  to  the  lower  boundary. 
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Figure  13.  Distribution-Function  Plots  for  Rain, 
Type  RjW  ,  for  Three  Different  Values  of  Liquid- 
Water-Content.  The  plots  are  similar  to  those 
for  ice-crystals  and  snow,  except  that  only  one 
abscissa  diameter  scale  is  shown  which  extends 
to  6  mm 


^Generally  true  for  'din  and  large-snow,  but  exceptions  can  occur  for  appreciable 
lower  diameter  truncation  in  the  case  of  small-snow  and  particularly  in  the  case 
of  ice  crystals. 
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8.  CONCLUDING  K  CM  ARKS 


This  report  has  described  how  values  of  liquid-water-content  in  the  precipita 
tion  size-range  of  the  hydrometeors  were  determined  for  the  SAMS  missile  tra¬ 
jectories  from  AFCRI.-acquired  radar  and  aircraft  data.  The  accuracy  of  the 
liquid- water-content  values  were  estimated  to  a  first  approximation,  and  infor¬ 
mation  was  supplied  about  the  size  spectra  of  the  hydrometeors  as  assessed  from 
theoretical  distribution  functions. 

It  should  be  emphasized  that  substantial  liquid-water-content  also  exists  in 
the  cloud-size  range  of  the  hydrometeors,  in  which  the  droplets  or  particles  are 
smaller  than  about  80  microns  diameter  (equivalent  melted  diameter,  in  the  case 
of  ice  hydrometeors).  The  RARF  radars  at  Wallops  Island  do  not  generally  de¬ 
tect  the  presence  of  these  cloud-size  droplets  or  particles;  consequently  their 
numbers,  s.zes,  and  liquid-water-content  must  be  determined  from  aircraft  mea 
surements.  Aircraft  data  and  size-distribution  information  concerning  this 
"cloud-size  portion"  of  the  total  hydrometeor  spectrum  will  be  included  in  later 
reports. 
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List  of  Symbols 


Arabic  Symbol* 

a  constant  of  the  fall  velocity  equation  for  hydrometeors 

A  horizontal  sampling  area  of  an  instrument  that  measures  rainfall 

or  rainfall  rate 

b  exponent  of  the  fall  velocity  equation  for  hydrometeors 

d  vertical  distance  downward  within  the  hydrometeor  transition  zones 

D  diameter  of  water  drops,  also  used  in  the  general  sense  to  include 

the  equivalent  melted  diameter  of  ice  hydrometeors 

De  equivalent  melted  diameter  of  ice  hydrometeors 

Dj  drop  diametir,  for  classified  data 

De.  equivalent  melted  diameter,  for  classified  data 

D0  median  volume  diameter  (or  median  equivalent  melted  diameter)  of 

the  drops  or  ice  particles  of  a  given  population 

D*  modal  diameter  (or  modal  equivalent  melted  diameter)  of  the  dis¬ 

tribution  of  liquid-water-content;  that  is,  the  diameter  (class) 
corresponding  to  the  peak  value  of  liquid-water-content 

D^.  modal  diameter  (or  modal  equivalent  melted  diameter)  of  the  dis¬ 

tribution  of  the  radar  reflectivity  factor;  that  is,  the  diameter 
(class)  corresponding  to  the  peak  value  of  the  radar  reflectivity 
factor 

dEM  uncertainty  of  the  values  of  the  exponent  of  the  M  vs  Z  equation 

dK^  uncertainty  of  the  values  of  the  constant  of  the  M  vs  Z  equation 

dM  uncertainty  of  the  values  of  liquid-water-eontent  (in  the  precipita¬ 
tion  size  range  of  the  drops  or  particles) 

dZ  uncertainty  of  the  radar-measured  values  of  the  radar  reflectivity 

rm  factor 


Preceding  page  blank 
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List  of  Symbols 


dZ 


s 


e 


E 


h 

H 

T 

I 

k 

^subscript 

K 

KM 

Kp 

l 


L 

M 


M 

subscript 

N 


Nt 

N(Dj) 

N(De.) 

Nf(Dj) 

Ng(Di) 

N(f(D>] 


uncertainty;  that  is,  scatter,  of  the  radar  reflectivity  values  (used 
in  regression  analyses)  about  the  regression  line 

base  of  natural  logarithms 

exponent  of  the  Z  vs  P  equation 

exponent  of  the  M  vs  Z  equation 

exponent  of  the  P  vs  Z  equation 

vertical  distance  upward  within  the  hydrcmeteor  transition  zones 
vertical  thickness  of  the  hydrometeor  transition  zones 
radar  integration  signal,  in  absolute  units  of  cm'1 
radar  integration  signal,  decibel  values 
constant  of  the  M  vs  P  equation 

constant  of  the  M  vs  P  equation  for  hydrometeor  region  and  type, 
or  other  conditions,  identified  by  the  particular  subscript 

constant  of  the  Z  vs  P  equation 

constant  of  the  M  vs  Z  equation 

constant  of  the  P  vs  Z  equation 

approximate,  average  physical  size  of  the  ice  hydrometeors;  that 
is,  of  the  snow  aggregates  or  ice  crystals 

length 

liquid-water-content,  of  precipitation-size  hydrometeors 

distributed  liquid-water-content— distributed  according  to  the 
diameter,  or  equivalent  melted  diameter,  of  the  hydrometeor 
drops  or  particles  — per  millimeter  bandwidth 

liquid-water-content  for  the  hydrometeor  region  and  type,  or 
other  conditions,  identified  by  the  particular  subscript 

number  concentration  of  the  hydrometeor  drops  or  particles  per 
m^  per  mm  bandwidth 

constant  of  the  non-truncated  exponential-distribution-function 
specifying  the  variation  of  the  number  concentration  of  the  drops 
or  particles  with  the  diameter,  or  equivalent  melted  diameter 

total  number  of  the  hydrometeor  drops  or  particles  of  a  given 
population,  per  m* 

number  of  drops  of  a  classified  diameter  size,  Dj ,  per  cubic  meter 

number  of  snow  or  ice  particles  of  a  classified  equivalent-melted- 
diameter,  De. ,  per  cubic  meter 

number  of  drops  of  a  classified  diameter  size,  Dj ,  that  fall  across 
a  given  horizontal  area  in  a  given  time— dimensionless 

number  of  drops  of  a  classified  diameter  size,  Dj ,  that  are  in¬ 
cluded  in  a  given  volumetric  sample  — dimensionless 

any  general  distribution  function  that  specifies  the  number  concen¬ 
tration  of  water  drops  as  a  function  of  their  diameter 
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List  of  Symbols 


N[f(De)l 


subscript 


Rs 

T 


any  general  distribution  function  that  specifies  the  number  concen¬ 
tration  o,  snow  particles  or  ice  crystals  as  a  function  of  their 
equivalent  melted  diameter 

precipitation  rate,  general 

precipitation  rate  for  the  hydrometeor  region  and  type,  or  other 
conditions,  identified  by  the  particular  subscript 

path  distance  along  the  SAMS  missile  trajectory 

time 


V 


V 

s 


A'/. 

Z 


Z 


D 


>M 


'IX 


J\V 


fall  velocity  of  the  drops  or  ice  particles 

fall  velocity  of  a  given  size  class  of  hydrometeors  having  the  mid¬ 
diameter  Dj 

sampling  volume 

absolute  altitude,  above  sea  level 

incremental  difference  of  altitude 

radar  reflectivity  factor,  general 

distributed  radar  reflectivity  factor— distributed  according  to  the 
diameter,  or  equivalent  melted  diameter,  of  the  hydrometeor 
drops  or  particles  — per  millimeter  bandwidth 

radar  reflectivity  factor  for  ice  hydrometeors,  as  measured  by 
radar 

radar  reflectivity  factor  for  ice  hydrometeors,  as  computed  from 
size-distribution  data 

radar  reflectivity  factor  for  ice  hydrometeors,  as  computed  from 
the  assumption  of  a  theoretical  distribution  function 

radar  reflectivity  factor  for  water  hydrometeors,  as  measured 
b  radar 


Z^y  radar  reflectivity  factor  for  water  hydrometeors,  as  computed 

'  from  size  distribution  data 

Z\y-p  radar  reflectivity  factor  for  water  hydrometeors,  as  computed 

from  the  assumption  of  a  theoretical  distribution  function 


(Irrck  S\iiiIm)Ik 


(V 

^subscript 

'^subscript 


Pint 


■y 

.subsc  ript 


constant  of  the  exponential  distribution  function  of  Marshall  and 
I ’aimer 

altitude  change  of  the  ecipitation  rate  within  the  hvdrometeor 
transition  zone  in<  ated  bv  the  particular  subscript 

altitude  change  of  the  constant  of  the  M  vs  P  equation  within  the 
hvdrometeor  transition  zone  indicated  by  the  particular  subscript 

gamma  function  of  the  particular  number  "n". 

altitude  change  of  the  exponent  of  the  M  vs  P  equation  within  the 
hvdrometeor  transition  zone  indicated  by  the  particular  subscript 


List  of  Symbols 


L 

£  subscript 


V 

V 
X 
A 


rr 


a 


Pr 


1 


exponent  of  the  M  vs  I’  equation 

exponent  of  the  M  vs  P  equation  for  the  hydrometeor  region  and 
type,  or  other  conditions,  identified  by  the  particular  subscript 

radar  volume  reflectivity,  ir  absolute  units  of  cm*' 

radar  volume  reflectivity,  decibel  values 

rad.-* r  v.  ;  ■'  -■>  gth 

expire'  i,.il  "slope  factor"  in  the  distribution  function  for  the 
number  concentration  of  the  drops  or  particles 

circle  circumference  divided  bv  diameter 

density  of  water 

elevation  angle  of  the  SAMS  missile  trajectory  above  the  horizon¬ 
tal  plane 

standard  deviation  of  pulse  or  video  integrated  received  power  for 
single  independent  samples 


Other  Symbol* 


when  used  to  "underline"  any  given  Arabic  or  Greek  symbol, 
represents  the  decibel  equivalent  value  of  the  parameter;  for 
example,  ,X.=  10  log  X  . 


it  U.  S.  GOVERNMENT  PRINTING  OFFICE:  19 7**- - 600- 0  1  J- -  10 


86 


